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INTRODUCTION  AND  RESEARCH  OBJECTIVES 


This  report  describes  work  carried  out  at  the  Gas  Turbine  and 
Plasma  Dynamics  Laboratory  at  M.I.T.,  as  part  of  a  multi-investigator 
effort  on  current  problems  in  turbomachinery  fluid  dynamics.  Support 
for  this  program  is  provided  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  Number  F49620-82-K-0002,  Dr.  J.D.  Wilson, 
Program  Manager. 

The  present  report  gives  a  summary  of  the  work  for  the  period 
10/1/81  -  9/30/82.  For  further  details  and  background,  the 
referenced  reports,  publications  and  the  previous  Final  Report 
(covering  the  period  10/1/79  -  9/30/81) 1  should  be  consulted. 

Within  the  general  topic,  four  separate  tasks  are  specified. 
These  have  been  described  in  detail  in  Reference  1,  but  they  are,  in 
brief ; 

1.  Investigation  of  fan  and  compressor  design  point  fluid 
dynamics 

2.  Studies  of  compressor  stability  enhancement 

3.  Fluid  mechanics  of  gas  turbine  engine  operation  in  inlet 
flow  distortion 

4.  Investigations  of  three-dimensional  flows  in  highly  loaded 
turbomachines 

In  addition  to  these  tasks,  the  multi-investigator  contract  also 
encompasses  the  Air  Force  Research  in  Aero  Propulsion  Technology 
program.  The  work  carried  out  in  each  of  the  tasks  will  be  described 
in  the  next  section. 

1.  E.M.  Greitzer  et  al.,  AFOSR  TR-82-0027,  Final  Report,  10/79-9/81, 
on  "Current  Problems  in  Turbomachinery  Fluid  Dynamics" 
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TASK  I:  INVESTIGATION  OF  FAN  AND  COMPRESSOR  DESIGN  POINT  FLUID  DYNAMICS 

TASK  IA:  INVERSE  DESIGN 


Task  Status 


During  the  last  two  years  a  procedure  which  uses  an  input  blade 
static  pressure  distribution  to  construct  a  geometry  for  supersonic 
cascades  has  been  produced.  This  procedure  assumed  inviscid  flow  in 
the  cascade.  Continued  progress,  for  either  two  or  three  dimensional 
design  schemes,  requires  that  a  definition  of  a  "better"  blade  row 
geometry  be  adopted,  and  we  chose  to  define  better  as  a  lower 
stagnation  pressure  loss  for  the  same  row  static  pressure  rise  .  The 
procedures  previously  developed  gave  us  the  ability  to  produce  a  blade 
row  design  for  a  specified  blade  static  pressure,  but  were  unworkable 
for  producing  better  cascade  geometries  for  two  reasons.  First,  for 
supersonic  cascades,  boundary  layer  growth  interacts  strongly  with  the 
inviscid  core  flow  and  makes  it  quite  difficult  to  design  realistic 
profiles  when  viscous  effects  are  not  considered.  Second,  for  low 
stagnation  pressure  loss  blade  designs,  the  numerical  errors  in 
conservation  of  stagnation  pressure  became  comparable  with  the 
physical  shock  loss. 

We  felt  that,  since  a  number  of  procedures  to  add  viscous  flow 
effects  to  an  inviscid  calculation  existed,  but  that  little  or  no 
recent  progress  had  been  made  toward  improving  conservation  properties 
of  time-marching  algorithms,  it  was  most  important  to  concentrate  on 
the  problem  of  stagnation  pressure  errors  in  numerical  algorithms. 

This  problem  turned  out  to  be  much  harder  to  understand  and  cure  than 
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had  been  anticipated,  and  nearly  all  our  efforts  this  year  have  gone 
into  solving  this  under-appreciated  problem. 

This  effort  has  largely  been  sucessful,  and  our  present 
algorithms  have  reduced  the  stagnation  pressure  error  from  an 
intolerable  10%  to  15%  to  about  1%.  At  the  present  time,  we  have  what 
we  consider  to  be  the  best  possible  second  order  accurate  numerical 
scheme  for  our  purposes,  and  we  are  using  it  to  design  shock-free 
supercritical  stators  in  order  to  demonstrate  that  we  have  solved  the 
stagnation  pressure  error  problem.  Our  next  step  will  be  to  combine 
this  improved  numerical  scheme  with  a  new  boundary  layer  code  under 
development  for  NASA  Lewis  so  that  optimization  calculations  including 
viscous-inviscid  interaction  may  be  done  for  the  design  task. 

Summary  of  Effort  for  the  Current  Reporting  Period 

The  major  effort  for  this  reporting  period  was  the  search  for 
sources  of  stagnation  pressure  errors  in  time-marching  calculations. 
Two  important  sources  of  error  were  identified  and  corrected  in  the 
current  algorithms.  The  first  error  source  was,  roughly  speaking, 
that  the  algorithm  solutions  were  sensitive  to  the  shape  of  the  finite 
difference  grid  in  a  way  distinct  from  the  fineness  of  the  grid.  In 
finite  element  terminology,  the  algorithms  did  not  maintain  geometric 
area  conservation.  The  second  error  source  was  due  to  the  details  of 
specifying  the  solid  wall  boundary  conditions. 

Correction  of  errors  due  to  grid  shape  is  an  important  step 
toward  solution  of  the  design  problem,  but  errors  due  to  specification 
of  solid  wall  boundaries  are  unimportant  for  the  projected 
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viscous-inviscid  interaction  calculation  since  the  flow  in  the  solid 
wall  region  will  be  calculated  by  a  boundary  layer  solver. 

The  search  for  these  error  sources  demonstrated  that  nearly  all 
popular  time-marching  schemes  have  identical  interior  point 
properties,  and  that  finite  difference  solutions  to  problems  posed  in 
strong  conservation  law  form  can  be  made  to  be  identical  to  so-called 
finite  volume  formulations.  This  work  was  collected  in  an  AIAA  paper, 
A1AA-83-0257  included  as  Appendix  A,  which  will  be  presented  in 
January  1983. 

Efforts  to  understand  these  stagnation  point  errors  did  not 
produce  a  dramatic  new  algorithm,  but  did  demonstrate  how  to  optimize 
existing  algorithms  to  produce  significantly  better  results.  The  most 
severe  test  of  our  new  numerical  procedures  is  the  shock-free 
supercritical  stator  problem  where  the  flow  solution  must  be  nearly 
isentropic  at  every  point  in  order  to  produce  shock-free  solutions. 

The  test  case  chosen  was  the  design  due  to  Sanz,  and  comparisons 
between  predicted  and  computed  surface  Mach  numbers  are  shown  in 
figures  12  and  14a  of  Appendix  A.  The  level  of  stagnation  pressure 
errors  is  shown  in  figure  13  and  14b.  Older  procedures  predicted  a 
shock  on  the  suction  surface  and  a  stagnation  pressure  error  greater 
than  15%.  The  present  results  are  extremely  good  for  an  Euler  flow 
calculation  matched  against  a  potential  flow  solver. 


In  order  to  demonstrate  that  it  is  feasible  to  optimize  cascade 


blade  shapes  using  the  procedures  developed,  three  projects  are 
planned  for  the  next  reporting  period. 

1 )  Design  a  supercritical  stator  similar  to  the  Sanz  blade 
with  a  smaller  gap  to  chord  ratio.  The  hodograph  procedure 
used  by  Sanz  often  breaks  down  as  the  blade  spacing  is 
reduced.  This  project  will  complete  the  thesis  work  of 
Tong. 

2)  Couple  the  present  numerical  code  to  the  viscous- inviscid 
interaction  procedure  under  development  for  NASA  Lewis  to 
incorporate  viscous  effects  into  design  procedure. 

3)  Design  an  optimum  two-dimensional  supersonic  compressor 
cascade  on  the  basis  of  minimum  stagnation  pressure  loss 
for  shock  system  loss  plus  wake  loss. 


I 
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APPENDIX  A 

SOLUTION  PROCEDURES  FOR  ACCURATE  NUMERICAL 
SIMULATIONS  OF  FLOW  IN  TURBOMACHINERT  CASCADES 

W.  T.  THOMPKINS  Jr.* 

S.  S.  TONG** 

R.  H.  BUSH** 

W.  J.  USAB  Jr.** 

R.  J.  G.  NORTON*** 

Massachusetts  Institute  of  Technology 
Cambridge.  Massachusetts 


Abstract 


INTRODUCTION 


For  several  years  the  authors  have  been 
involved  in  evaluating  and  developing  numerical 
simulation  schemes  for  compressible,  two 
dimensional  inviscid  or  viscous  flows  in 
turbomachinery  cascades.  Numerical  schemes 
considered  were  all  originally  classified  as 
time-marching  schemes  and  include:  1)  implicit 
approximate  factorization  schemes  like  those  of 
Beam  and  Warming,  2)  explicit  schemes  due  to 
MacCormack,  1)  explicit  central  difference  schesws 
due  to  Jameson,  Schmidt  and  Turkal  and  Rizzi,  and 
4)  the  multi-grid  scheme  of  Mi.  As  we  have 
developed  these  schemes  we  have  cosw  to  believe 
that  the  accuracy  of  computational  results  is 
relatively  insensitive  to  the  numerical  algorithm 
chosen  but  highly  sensitive  to  implementation 
details  such  as  boundary  conditions,  consistent 
flux  balancing,  grid  resolution  and  numerical 
smoothing.  In  order  to  illustrate  our  viewpoint, 
we  present  an  examination  of  the  relationship 
between  a  flux  balancing  interpretation  of  the 
control  volume  conservation  laws  and  various  finite 
difference  formulations  and  comparisons  of  the 
performance  of  these  schemes  on  three  test 
problems:  Nl's  bump  in  a  channel,  a  supersonic 
nozzle,  and  flow  in  a  supercritical  compressor 
cascade.  As  a  final  teat  example,  a  comparison  is 
made  between  measured  and  computed  flow,  using  a 
compressible,  Navier-Stokes  solver,  for  a  high 
speed  turbine  cascade. 


Finite  Difference  Operator  Definitions 
"j.*  *  f<ttJ*l.k  *  “j-1,k> 

A  uj,k  *  |('*J*1,k  *  »j.k» 

ui.k  ‘  *  uj-1,k> 

*3  uj,k  *  u}*i,k  *  uj.k 

A  aJ»k  *  u3.*  *  ui-».k 
s}/2  uJ.k  •  “3-4, k-  “3-i.k 


Note  that  all  coordinate  positions  at  non-integer 
mesh  spacing  are  to  be  defined  by  simple  averages 
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For  several  years  the  authors  have  been 
involved  in  evaluating  and  developing  numerical 
methods  for  compressible,  inviscid  or  viscous, 
internal  flow  problems.  The  methods  discussed  in 
this  paper  include  ( l )  implicit  approximation 
factorization  schemes  due  to  Beam  and  Warming  [1], 
(2)  explicit  schemas  due  to  MacCormack  [2],  (3)  an 
explicit  central  difference  scheme  suggested  by 
Jameson,  Schmidt  and  Turkel  [3]  and  Rizzi  [4),  and 
(4)  the  multi-grid  schema  due  to  Ni  (S).  Most  of 
these  methods  were  originally  formulated  as 
time -marching  methods  although  all  are  now  run  as 
pseudo-time-siarching  methods. 


Throughout  these  evaluations  we  have  been 
interested  primarily  in  the  accuracy  of  the  steady 
state  solutions  rather  than  the  coagmter  run  time 
or  number  of  iterations  to  convergence.  As  we  have 
developed  these  methods  we  have  cm  to  believe 
that  the  accuracy  of  the  computational  results  is 
relatively  insensitive  to  the  numerical  schema 
chosen,  but  it  is  highly  sensitive  to  implemen¬ 
tation  details  such  as  boundary  conditions, 
consistent  flux  balancing,  grid  resolution  and 
numerical  smoothing. 


In  order  to  illustrate  our  viewpoint  we  first 
present  an  examination  of  the  relationship  between 
a  flux  balancing  or  finite  volume  Interpretation  of 
the  control  volume  conservation  laws  and  a  finite 
difference  formulation  in  strong  conservation  law 
form,  while  the  development  hare  is  perhaps  not 
unique  or  entirely  new,  it  is  essential  for  under¬ 
standing  the  relationships  between  various 
numerical  schemes.  Next  we  present  a  development, 
from  a  flux  balance  viewpoint,  of  the  several 
schemes  which  emphasizes  their  basic  similarity 
rather  than  their  computational  differences.  The 
conclusion  hare  is  that  all  these  schemes  should 
produce  nearly  the  sai  s  steady  state  solutions  when 
the  seme  boundary  conditions  and  smoothing 
operators  are  applied.  Results  from  three  inviscid 
computational  test  problems  are  then  presented  in 
order  to  demonstrate  the  validity  of  this 
conclusion.  Lastly,  results  from  computations 
using  a  full  Navier-Stokes  version  of  the 
Beam-Warming  algorithm  are  compared  to  experimental 
results  for  a  high  speed,  high  turning  turbine 
cascade.  These  results  illustrate  that  a  consis¬ 
tent  application  of  the  flux  balancing  analysis  can 
produce  accurate  solutions  even  in  extrema 
geometries  with  rsthar  modest  cost.  The  number  of 
iterations  to  full  convergence  is  about  500  for  a 
grid  containing  100>50  nodes. 
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Salatlonshlp  Between  finite  Difference  and 


Consider  the  two-dimensional  Eular  aquation 
axpraaaad  in  conssrvation  law  form 


—  *  2£  +  i£  «  o 

it  ix  Jy 


ouv 

ou(£+P/o ) 


PV 

OUV 

0f3  *  p 

ov(S+P/o ) 


and  a  coordinata  system  transformation. 


C  •  C(x,y>  n  ■  «(x,y)  t  »  t  .  (2) 

Tha  strong  conservation  law  fora  undar  this 
transformation  is  ratainad  as,  rafaranca  (61 : 

_ifs^  *  _»(e.  ii  ♦  a  !i  1  ♦  sa  ♦  «  !a )  .  „ 

7t\jJ  ie^j  ix  J  iy  J  in^J  ix  J  ly  J  (3) 
where  J  is  tha  transformation  Jacobian 


In  shorthand  notation,  aquation  (1)  becomes i 


°.t  *  r.(  *  °.n  “  0 


n  /J  , 

,y 


c  y/j  . 


'  \x/J  ' 


«,*/J  ' 


and  the  inverse  Jacobian  is 

J“1  *  X,C  y,n  ‘  x,n  y.C  *  VJ  • 

The  Jacobians  express  the  ratio  of  call  areas 
between  tha  different  representations. 

The  transformation  metrics  (C,x  E,y 
n(y)  are  most  conveniently  obtained, 
suggested  by  stager  [J] ,  from  a  finii  ..  ference 
evaluation  of  those  derivatives  in  o  tational 
space  using  equations  (6).  This  proi  uaes  only 
tha  fact  that  the  coordinate  transfox  .Jn  is  a 
one-to-one  napping.  If  we  consider  i  iese 
derivatives  to  be  evaluated  by  cents'  iuite 
differences,  we  have,  assuming  AC  -  . 

(Vjj.k  '  **  yj.*  *  (7) 

(Vjj.X  *  ~  xj,X  *  — ~  (8) 

y^-^^vrj’u-  ,9) 

(Vji.x  •  *3  *3.*  • 


•  *2  x. 


3  x3.* 


.  tK*M.k  ~  *1-1. K1 


J  ■  *.x  ",r  "  c.r  ",x 

Formal  relationships  exist  between  tha 
transformation  from  (x,y>  coordinates  to  (C,n) 
coordinates  and  tha  inverse  transformation  from 
( C.n)  coordinates  to  (x,y )  coordinates  which  are i 


If  we  consider  the  computational  cells  to  be 
polygons,  the  exact  area  of  a  cell  can  be 
determined,  see  [B] ,  from 


Area  -  £  xt  yw,  -  £  xle1  yi 

i-1  i-1 


where  N  is  tha  number  of  edge  nodes,  and 


i*l  if  i*1  i  N 

1  if  i*1  >  M 


Figure  1.  Relationship  between  transformation 
from  physical  to  computational  space 
and  computational  to  physical  space. 


Figure  :.  lode  connection  diagram  for  general 
computational  cell  polygon. 
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Mow  consider  a  possible  finite  difference 
representation  of  equation  { 5 ) , 


4j  •:.* 


-  i 


o  ; 


'j.k 


( 12) 

[  r 

on  a  standard  9  point  cell. 


j-i, it*  i 

j.kt-1 

j+l.k-H 

j-1, A 

j.k 

j+1,k 

j-1 ,1c- 1 

j.k-1 

Ji-i.k-1 

or 


‘j-1 ,k-l il 


( 14) 


This  form  of  the  continuity  equation  is  identical 
to  that  which  would  be  obtained  by  ac  integration 
of  the  control  Toluna  fluid  equations  around  the 
cell  area,  assuming  no  variation  of  (ou)  or  (ov) 
across  a  cell  face. 


figure  1.  Finite  difference  node  numbering  scheme. 


Again  using  only  the  continuity  equation  for 
clarity. 


Sow  replacing  the  metric  derivatives  with  their 
finite  difference  expressions,  equations  (7) 
through  (11),  and  writing  out  only  the  continuity 
equation  for  clarity,  we  have 


r  < 


Vj+1 ,k+l  ■  y j*l  ,k 


'  (0V)  •  xj*1.k-l)j 

-  -  yj-i,k-i) 

■  HS-mK’.-i  *  ’e^-’ 

e 

*  [*  W3,h*’(y^1’k<-1  ‘  y3-’ .*•») 

■J 


a 

TE 


Jadxdy  -  -  J 
area  - 


o(u  i*  *  v  iy)*n  dl 
control  surface 


(16) 


where  ix,  iy  are  unit  vectors  in  original 
Cartesian  coordinates, 

n  is  an  outward  normal  vector  to 

the  cell  face,  and 

di  is  length  along  the  face. 

Equation  (16),  when  expanded,  is  identical  to 
equation  (14),  and  shows  that  the  proper 
interpretation  of  a  metric  derivative  is  that  the 
ratio  of  that  derivative  to  the  Jacobian  is  the 
projected  area  of  the  face.  Thus  the  solution  of  a 
finite  difference  problem  in  strong  conservation 
law  form  in  a  transformed  space  is  equivalent  to 
the  solution  of  a  flux  sumation  algorithm 
expressed  in  physical  space. 


The  flux  sun  of  equation  (13)  could  be 
improved  by  electing  a  more  accurate 
approximation  for  the  integral  of  (ou)  or  (ov) 
over  a  face  than  simply  its  value  at  the  center 
of  the  face.  For  example,  a  piecewise  continuous 
interpretation  of  (ou)  is  equivalent  to 

j"(ou  ‘-„>*n  d(  -  (ii£(ou>j,)cj  yj#< 

*  (“k<®u>j,i()  4k  yj,k  •  '17> 

An  important  property  of  equation  ( 14 )  is  that 
a  uniform  solution  t(eu)j(j.  ,  (ov)^  ^  constant 
over  space)  remains  identically  uniform.  This 
property  is  generally  referred  to  as  area  conserva¬ 
tion  in  finite  element  analysis  and  its  importance 
to  maintaining  constant  freestream  solutions  was 
first  pointed  out  by  Steger  ;7).  a  recent  paper  by 
Hindman  ’9;  has  also  emphasised  this  property. 
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COMPARISON  OF  FLCX  BALANCE  FORM  0?  SEVERAL 

TIME-MARCHING  ALGORITHMS  FOR  INVT5CID  FLOW  The  flux  balance  optrttot  is: 


Csntsrsd  Dlffsrancs  Algorithms 

METHOD  A.  Implicit  Time-Msrchinq  Algorithm  of 
Bsan  «nd  Warming  tils 

The  first:  method  to  ba  considered  is  tha 
-Ur a  starching  algorithm  ineroducad  by  Beam  and 
Waning  which  is  unconditionally  ttabla,  in  tha 
linaar  problem,  for  two  dimensional  flows: 


“  “  jc  ■  Ljwav  Gj,k)  <  <18> 

share  l{  and  t.Jj  are  finite  diffaranca 
operators  whose  fora  datarainas  tha  order  of 
accuracy  of  the  time-marching  algorithm  and  L,BwkV 
is  a  finite  diffaranca  smoothing  operator. 

Solution  Procedure 


RC"j,k  ■  «“/2(FCH)",lc  ♦  SS/2<CCH)5,<  .  («) 

where 

(rCHrji*,k  “  (ui  r".k)(“j  *k/2  yj.k) 

-  [hj  Sj.kj^j  <k/2  Xj.x| 

(GCH)"(Wj  -  -  ^uk  »j/2  yj.k) 

*  (“it  Gj,k)  (»k  s°/2  *J,k)  * 

and  the  algorithm  bacootas 

(I  *  4tj,ktC,(I  *  ^j.k^n^-jA  "  ®j.k] 

-  -  Atj,k((RCH)5,k  -  Lawkv  Oja)  •  (23> 


Step  1 : 


(l  *  At  Lc](f0j,k)  *  -  »*(»«},„  -  WWAvSj.k)  (19) 
Step  2: 

fl  ♦  At  L$j  A0"fk  -  Aojfk 


(20) 


Step  3: 


~nv1  ~n  _  ~n 
cJ,k  *  u3,k  “i.k 


(21) 


METHOD  B.  Modified  Beam  and  Waraino  Algorithm 

In  order  to  improve  tha  accuracy  of  the  flux 
balance  operator  in  skewed  mash  systems,  a  variety 
of  modified  Beam  and  Harming  algorithms  have  bean 
developed.  The  simplest  modified  operator,  which 
corresponds  to  several  other  algorithms,  is  to 
balance  the  flux  on  tha  call  shown  below. 


j-1.k*1 

j,k*l 

3*1 ,k*1 

j“t.k*i 

1 

I 

mm 

j-t,k 

_ _ _ J 

:v!;lyXv!i 

H 

3-t.k-i 

■ 

}*1.k 

A 

1 

•n 

I 

1 

j.k-1 

li*i,k-l 

1 

Figure  4.  Half-spacing  mash  cell  in  computational 

space. 


METHOD  C.  Explicit  Tima  Integration  Scheme  of 
Jameson  (31  and  Rlrrl  (4)  : 

Tha  Runga-Eutta  time  stepping  approach  can  ba 

expressed  as: 


-  -  RCH^k  *  l^v  o$,k  -  -  po^,k  (24> 


Step 

_1: 

S(0)  , 

5l,k 

(25) 

Step 

_2: 

-(1) 

0  » 

5<0)- 

(26) 

Step 

_3: 

5(2>  . 

5'”-  j 

(%, 

(27) 

Step 

_4: 

*13) 

0  » 

S'”  . 

M2) 

^  1tk  ™ 

(2B) 

Step 

J5: 

S'4’  . 

S'” 

- 

fit) 

^PD,0>*  2PU:1>-  2J>S(2>-  po(3>) 

(29) 

stea. 

_6 

(30) 
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Thi>  method  is  usually  run  with  the 
smoothing  operator  in  aach  stap  formad  as 

-n  » ( i ) 

l»RKAV  u  rathar  than  Irxav  o  • 

If  staps  2,  3.  4,  S  aach  converge  independent- 
ly,  than  the  staady  stata  solution  is  identical  to 
that  of  tha  modified  Beam  and  Wanting,  method  B,  if 
LBWAV  *  lrxav>  and  tha  solution  is  independent 
of  4tjfk. 

METHOD  p.  Flux  Balance  Form  of  MacCornack's 
Method  [21: 

A  flux  balance  font  of  MacContack's  method 
may  be  expressed  as: 


Step  1,  Predictor: 


“j,k  ”  “j,k  “  4tj,k  H^j+l.kW-yl  ”  Gj-M,kpl'xlj 
"  (Fj.kKTr3  -  Gj,kK*x3  j 

♦  (“  *  Gj,k+1PLx4j 

-  (-  Fj,kPty2  +  Gj.kK-xa  ) 

-  r,n 

*  Lmav  ^j.k 


where  the  projected  lengths  ara 
PLy,  -  uj  4k/2  yj<k  -  Yj^.k+i  -  Yj+*,k-fc 


PLx4  -  uk  tV2  xj(k  ’  xjH.*H  _  *1-i.k-4 


(31) 


Stap  3  nay  be  rewritten  as 

*n*  1  ;n 

Uj,k  ’  -j,k 

*  (<rj*1,k  *  Fj,k><“j  sk/2  Yj,k> 

"  {Gj+1,k  *  Gj,k){llj  4k/2  xj,kl] 

"  |<pj,k  *  fj-l.kXkj  *k/2  Yj,k> 

-  t5j,k  *  ®j-t.k,t“j  4k/2  xJ,kl) 

*  I"  {pj,k*1  *  fj.kHi'k  4  j/2  yj.kl 

*  (GJ,k*1  *  Gj,kJf“k  4  j/2  xj,kJ] 

“  (-  <pj,k  ♦  pj,k-lH»k  4 j/2  *j,k> 

*  (Gj,k  *  Gj,k— 1><>*k  4 j/2  xj,k^j 

-  U^v  {Oj.k  *  “j.k1 


(34) 


This  font  of  MacContack's  method  is  area 
preserving  for  aach  step  (31  or  32)  and  offers 
the  possibility  of  steady  state  solutions  which 
are  independent  of  it.  When  the  predictor  and 
corrector  step  both  converge,  which  we  generally 
find  to  be  the  case,  then  the  staady  state  form 
of  equation  { 34 )  is  t 


«CHj.k  -  I-WV  “j.k  -  0  ■ 


(35) 


as  with  tha  modified  Beam  and  Warming  and  fourth 
order  Runge-Kutta  schemes,  methods  B  and  C. 


Stap  2,  Corrector: 


Comments  on  Cantered  Difference  Algorithms 


“j!k  *  “j,k  "  4t j ,k |  (^j.k^yl  *  Gj,k«-xl] 

L'  1 

~  (rj-1,kPiy3  *  Gj-1,kPLx3  ) 

*  I*  Fj,kPLy4  *  Gj,k«nc4) 

*  I"  Pj.k-I^J  *  G j , k- 1  ^x2  ] 

*  ^MAV  “j.kl  *  <32) 


For  each  of  tha  methods  considered.  A,  B,  c  and 
0,  the  staady  state  solution  desired  is  that  the 
cantered  differenced  flux  balancing  residual 
(RC)S/k  or  ,RCH*?,k  z**°‘  For  »H 
these  schemes  at  interior  points,  this  criterion  is 
possible  only  when  the  artificial  smoothing 
operators  are  neglected.  The  order  of  accuracy  of 
all  these  methods  is  the  same  and  aach  should  be 
expected  to  produce  nearly  the  same  steady  state 
solution  except  for  differences  in  smoothing 
requirements  or  boundary  conditions.  As  we  shall 
see  in  the  test  examples ,  this  conclusion  can  be 
demonstrated  by  numerical  experiments. 


Mon-Centered  Difference  Algorithms 


“n-M  l  f-n  *#•  ) 

Jj,k  •  1  l°j,k  *  °j,kl 


METHOD  S.  MacCornack's  Method  on  a  Mon-Centered 
Cell  Basis: 


A  version  of  NacCormack's  method  which  does  a 
flux  sua  on  the  cell  structure  shown  in  Figure  5 
has  been  developed  by  Tong  MO).  This  form  is  a 
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two  or  three  dimensional  extension  zi  the  strong 
conservation  lew  vsrs.  cn  of  .MacCcmac'. •  s  method 
proposed  by  Hi.ndaan  '  9 ' .  This  senses  ises  the  same 
basis  nine  point  call  as  m  Figure  ;•  ;r  4,  but  it 
coaputas  a  flux  sun  ovit  each  of  t.ta  four  interior 
cells. 


nj/<  *  °j,k  ~  4tj,kf<w<c)  ~  Vacav  °j,k] 

‘  7  (3 

Step  2; 

°3%  “  °$,k  -  4tj,k[ (WIC)  -  Vav  05 #k] 

1  (3 


Step 

_3: 

„!»♦  1 

1  I  1 

Jj,k  - 

2  !° 

„n+1 

n 

°j,k  - 

°j.k 

♦  <RMC>;_4,k_4  -  WexvtOj.S  ♦  -Jj.k-' 

where 

mri”  .  «1  +  «k  ««*:•  = 


«i  ♦  «k  «ar:-:.k 


,FM  •".k  “  rj.kskyj.k  ~  °J,k4k*j,v.  . 


<G"~>".k  *  _  rj.k4jyj.k  *  Gi,k4jx;,k 


hjt1/2,kt1/2  !L*  the  area  of  4-point  cell. 

This  torn  of  NacCoreack's  method  also  obeys 
geometric  area  conservation  on  the  i-ooint  cells 
with  the  possibility  of  a  steady  state  solution 
that  is  Independent  of  at.  Both  the  predictor 
and  corrector  steps  can  individually  converge. 
The  predictor  and  corrector  steps  express  flux 
conservation  on  the  4-point  cells  and,  if  each 
(BMC) j>k  converges  to  zero,  then  we  have  global 
flux  conservation. 

METHOD  f.  Mi's  Lax-Wendroff  Method  '!'  « 

Vi's  method  may  be  expressed  it  a  one  step 
fore  as: 


In  ,00  _  t 

°3.k  ’  aj,k  *  **J,k‘»;/2  “k/2 


‘  —MIAV  3Jj,k 


(««ll±4.kt4 

x  / 

’  *7**!**'  ^(Hl.ki*  -  «w(axj5**.k  <  <“> 

(n*lj,k±* 

•  (“k  fj,k)(«k  Fj,k)  -  (uk  Sj,k)  (*k  xj,k]  «  '4S» 

(GOT)"t*'k 

‘  ■  (“j  Fj,k](4j  yj.k]  *  (“j  G",k)(‘j  xj,k)  <«> 

This  representation  of  Hi  schena,  which  is 
not  the  coeputational  fore,  is  a  Lax-Wendroff 
type  aethod  which  does  a  flux  sub  on  4-point 
cells.  A  possible  solution  to  the  steady  state 

form  of  aquation  (43)  is  that  each  (BHD j±i/2,)c±i/2 

becomes  zero  in  the  absence  of  artificial 
smoothing. 

Comments  on  Hon-Centered  Schemes 

Both  the  Tong  and  the  Hi  non-centered  schemes 
compute  flux  balances  on  the  4-point  cells  and 
each  solution  can  be  consistent  with  these 

balances,  (EMC) j41/2,k±1/2  or  < J*1/2,ki 1/2  ' 

becoming  zero  in  the  absence  of  smoothing.  Thus  the 
Tong  and  Ni  schemes  should  produce  approximately 
the  same  steady-state  solution,  if  the  same 
boundary  conditions  are  applied.  The  two  steady 
state  flux  balance  operators  are  not  Identical.  He 
would  expect  Nl's  scheme  to  produce  a  more  accurate 
result  in  the  absence  of  artificial  smoothing,  but 
his  method  nay  require  more  smoothing  than  Tong's 
method . 

nmsCID  TEST  EXAMPLES 

Three  configurations  were  chosen  as  test 
examples .  The  first  case  was  the  bump  in  a  channel 
problem  introduced  by  Ni  [51  j  the  second  was  a 
supersonic  nozzle  with  a  rapid  isentropic 
expansion:  and  the  third  was  a  supercritical  stator 
designed  by  Sanz  [ill.  For  each  problem,  an  exact 
steady  state  solution  would  be  isentropic  with 
constant  stagnation  pressure,  stagnation 
temperature  and  mass  flow  rate.  Methods  are 
compared  on  the  basis  of  predicted  Mach  number 
distributions  and  stagnation  pressure  loss  over  the 
domain,  ror  all  methods  on  all  problems,  the  mass 
flow  rate  was  conserved  to  within  0.3%. 

Boundary  Conditions 

Inflow/outflow  boundary  conditions  were  of  the 
non-reflecting  characteristic  or  extrapolation 
type,  but  implementation  details  differed  greatly 
from  method  to  method.  We  believe  that  the  results 
presented  are  as  free  as  possible  from  contamina¬ 
tion  due  to  these  conditions. 
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icsvef.gescs  criterion 


It  was  impossible  to  maintain  axactly  the  sama 
hard  wall  boundary  eondltiona  batwaan  tha  schemes. 
Methods  A  and  B,  which  wara  davalopad  atrictly  lor 
viscous  flow  by  tha  authors,  used  only  a  simple 
pressure  and  parallal  velocity  extrapolation  for 
the  inviscid  case.  These  conditions  wara  aufficiant 
for  tha  first  two  cases  but  not  the  supercritical 
stator  casa.  The  proper  conditions  for  Ni's  method 
are  not  clear,  Tha  first  two  examples  for  Ni's 
aethod  predicted  surface  quantities  fro*  tha 
interior  point  schaaa  and  than  corrected  thasa 
quantities  by  a  siapla  one— dimensional  wave  to 
maintain  flow  tangency. 

Methods  C,  0,  and  C  used  boundary  procedures 
which  consistently  satisfied 

normal  Momentum 

i£  -  £S_  (47) 

In  R 

L  I  1  * 

whara  q  *  u  ♦  v 

and  R  is  tha  radius  of  curvature  of  the  surface 
streaaviae  onaentua 

>£  >  -  0  (48) 

la  Is 

where  s  is  the  arc  length  along  surface 
flow  tan gen cy 

v/u«  4jy  /  4jx  (49) 

at  each  point  on  the  boundary.  To  satislfy  all 
three  conditions  at  each  point  on  the  boundary 
requires  an  iterative  evaluation  of  the  boundary 
pressures  at  each  time-march ing  step.  This 
procedure  is  costly  but  is  essential  for 
simulations  with  minimum  stagnation  pressure  error. 


artificial  Smoothing  Operators 

For  methods  a.  8,  D,  and  F,  the  smoothing 
operators  ware  generally  as  proposed  by  their 


authors  with  the  smoothing  coefficients 
for  each  calculation.  For  methods  C,  0, 
calculation  procedure  wasi 

"tunsd-up" 
and  E,  tha 

U*  -  ;i  -  R8FB)Un 

(50) 

o**  »  u  *  sx;u* 

(51) 

0n»l  .  ;i  *  s„)U** 

(52) 

where  RSFB  is  the  flux  balance  operator 
<!x  ,  Sy  are  the  smoothing  operators 

end 

SxU  •  ft  Uj*,*  Uj.^Ujl 

(53) 

f-  f’  ft*  1 

1*f2  l  4x/4x*tnj 

(54) 

f,  -  0.1  fj  -  200. 

For  methods  C  and  0,  these  smoothing  operators  were 
applied  after  each  step,  but  for  aethod  X  they  were 
applied  only  after  the  last  step. 

The  steady  state  solutions  are  not  independent 
of  the  smoothing  operators,  but  every  effort  was 
made  to  raduee  their  effect. 


Tzz  all  --  the  methods  except  Tong’s,  tha 
convergence  criterion  was  easy  to  establish.  The 
computed  solutions  were  run  until  the  steady  state 
flux  balance  operator  plus  the  smoothing  operator 
converged  to  a  specified  value,  1.0  X  10“* 
maximum  at  any  point  and  1.0  x  10*6  root  mean 
square  value. 

Is  Tong's  method,  as  with  any  MacCormack  type 
method,  the  order  of  differencing  in  predictor  and 
corrector  steps  is  not  unique.  The  double  sawtooth 
solutions  s .tenon  to  central  difference  schemas  can 
be  reduced  by  cyclically  permuting  the  operator 
sequence.  The  disadvantage  of  this  process  plus  the 
smoothing  sequence  adopted  is  that  we  lose  the 
ability  to  drive  a  steady  state  residual  below 
about  1.3  x  1 0"4  maximum  local  value  and  about 
1.0  x  10-5  root  mean  square  value.  It  is  thus 
difficult  to  determine  when  to  stop  iterating  on  a 
solution,  and  we  must  supply  some  heuristic 
principle  to  determine  when  to  stop.  The  advantage 
of  this  process  is  that  the  stagnation  pressure 
errors  in  the  "converged"  solution  are  greatly 
reduced. 


Si  Sumo  Test  Case,  Subsonic  Flow 

411  of  the  methods  produced  acceptably  accurate 
solutions  for  the  channel  bump  problem  especially 
in  terns  at  symmetry  and  Mach  number  predictions. 
The  grid  used  contained  65X17  points  and  is 
illustrated  in  Figure  S  along  with  a  typical  Mach 
number  contour  plot.  The  predicted  top  and  bottom 
wall  Mach  numbers  for  each  method  are  shown  in 
Figure  6.  The  stagnation  pressure  loss  for  each  of 
the  methods  is  shown  in  Figure  7,  \11  of  the 

centered  schemes  produced  essentially  the  same 
pressure  loss  of  about  2%  max  4Pt/Pt.  The 
non-ceatered  schemes  produced  essentially  no  loss. 
On  balance,  ell  the  schemes  performed  adequately  on 
tha  busp  case. 


Supersonic  Noxrle  Test  Case 

The  supersonic  noxzle  test  case  geometry 
and  grid  is  shown  in  Figure  8.  This  grid  has  65X33 
points.  The  inflow  Mach  number  is  2.0,  and  an 
idealized  wave  diagram  is  also  shown  in  Figure  8. 
The  wave  diagram  shows  that  the  wall  outflow  Mach 
number  should  be  2.83,  and  that  tha  Mach  number 
should  be  uniform  in  region  1 .  The  idealized  weve 
diagram  does  r.ot  account  for  curvature  of  tha 
character:.* tics  as  the  wave  system  from  the  lower 
wall  inters— s  with  the  upper  wall  wave  system. 
Under  this  assumption  tha  axit  cantarlins  Mach 
number  should  be  2.625.  The  results  for  aach  of 
tha  schemes  is  shown  in  Flgurs  9  in  terms  of  Mach 
nua&ar  centsurs.  The  correct  weve  system  is 
adequate.-/  predicted  by  all  schemes.  Ths 
stagnation  treasure  loss  contours  for  this  tast 
casa  srs  all  shown  in  Figure  10.  The  centered 
schemes  all  have  unacceptable  losses  In  tha  range 
of  34  to  44.  Ni's  scheme  has  a  small  loss  whils 
ths  tor.-car.tsred  MacCormack  has  virtually  no  loss. 
It  is  likely  that  a  bettsr  hard  wall  boundary 
condition  treatment  with  1 ' s  method  would 
eliminate  ths  loss  shown. 


t 
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Tha  ijrwat  bttMU  eh*  4  centered 
schemes  ia  quit*  r— rkabla;  In  fact,  a*  would  b* 
expected  from  eh*  analysis,  all  schaees  produo* 
nearly  eh*  lua  staady  seat*  solution*  •  Th* 
non-c*ne*r*d  MacCormack  and  Hi's  schaaa  produe* 
quit*  slnllar  r* suits. 

Supercritical  Stator  Tsst  Cas* 

A  meh  nor*  dasundiag  geometry  1*  th* 
supercritical  stator  design  pr***nt*d  by  Sans  (11], 
Th*  grid  for  this  test  cas*  is  shown  In  Figure  11. 

A  sispl*  sheared  grid  was  used  for  simplicity  which 
required  a  fin*  grid  near  th*  blad*  leading  *dg* 
for  accurate  simulations. 

Th*  surface  Mach  number  predictions  tor 
th*  supercritical  stator  using  th*  Tong1*  non- 
c entered  MacCoraack  type  method  are  shown  in  Figure 
12.  Th*  agreement  is  satisfactory  except  at  8% 
chord  on  the  suction  surface.  Slight  variations  in 
th*  specified  downstream  pressure  specified 
produced  overshoots  or  undershoots  in  this  region. 
The  saagile  calculation  is  for  th*  downstream 
pressure  specified  by  Sans.  Me  are  simply  unsure 
of  the  reason  for  this  discrepancy.  The  stagnation 
pressure  loss  for  this  calculation  is  shown  ia 
Figure  13.  An  unacceptably  high  lose  of  148  exists 
locally  at  th*  stagnation  point  and  a  maximum  loss 
of  4*  exists  on  the  suction  surface. 

The  centered  MacCormack  scheme  performed  nearly 
th*  aame  on  th*  Sans  blade  as  is  shown  in  Figure 
14.  Th*  stagnation  pressure  loss  at  the  leading 
edge  is  locally  14*  and  a  maxi  mum  of  48  in  th* 
remainder  of  th*  domain. 


VISCOUS  TURBIN*  CALCULATION  example 

While  th*  results  from  the  supercritical 
stator  geometry  illustrate  that  accurate  ruler 
equation  solutions  can  be  obtained,  the  situation 
is  less  clear  for  viscous  flow  problems.  To 
illustrate  our  level  of  computational  abilities  for 
these  flows,  a  calculation  of  th*  flew  ia  a  high 
speed,  high  tuning  turbine  cascade  is  presented. 
Th*  mmwrical  method  used  is  similar  to  th* 
approximate  factorisation  method  of  8 earn- Warming 
(1]  except  that  th*  flux  balance  operator  usas 
trapesoldal  rule  integration,  as  la  equation  (17), 
on  the  computational  cell  of  Figure  4.  Local 
values  of  At  are  used  to  produce  a  constant  local 
CFL  number  of  about  4.  Details  of  this  algorithm 
will  be  available  in  reference  (121.  solutions  for 
turbine  geometries  generally  require  about  500 
iterations  to  converge  to  machine  accuracy,  « 
decimal  digits. 

Th*  experimental  results  and  geometric  details 
for  th*  cascade  were  made  available  ia  reference 
(13],  and  th*  geometry  is  shown  in  Figure  IS.  flu 
turning  la  12*  degree*  with  an  outflow  Mach 
number  of  0.7S.  The  design  neynolds  number  ia 
300000.  All  calculations  assumed  ' — <  — -  flow. 

The  computational  results  are  presented  la 
Figures  1*  and  17.  The  first  figure  ooegaree 
predicted  blade  surface  preesur*  to  the 
e^arlmsntal  valuaa  for  a  laminar  aeleelation  at 
design  Reynolds  number.  Th*  grid  chosen  has  100x28 
points  aad  ia  nearly  orrhogneel  at  each  nods  point. 


Predicted  results  ere  quite  sensitive  to  th*  grid 
resolution.  Reasonably  accurate  surface  pressures 
were  predicted  everywhere  including  the  stagnation 
point  and  the  trailing  adga.  The  naxlrau*  stagnation 
pressure  error  outside  the  boundary  layer  was  about 
18.  Wa  are  as  yet  unable  to  explain  th*  discrepancy 
at  708  chord  on  tha  auction  surface  although  w* 
believe  it  to  be  linked  to  th*  trailing  edge  flow 
field. 

The  predicted  trailing  edge  flow  field  is 
Illustrated  in  Figure  17  in  terms  of  velocity 
direction  vectors.  Rile  calculation  shows  two  small 
Reparation  zones  comparable  in  alsa  to  th*  trailing 
edge  radius.  While  the  predicted  flow  field  at  the 
trailing  edge  seams  reasonable,  we  have  no  way  to 
verify  its  accuracy. 


CONCLUSIONS  AMD  DISCUSSION 

Tbs  presentation  of  tha  flea-marching 
schemes  ia  terms  of  a  flux  balancing  interpretation 
and  the  in  viscid  computational  asaaples  shows  that 
the  six  schemas  analysed  ex*  quite  similar  la  their 
steady  state  solution  properties.  While  th* 
computational  details  of  the  schemes  are  very 
different,  the  steady  state  eolations  of  th* 
central  difference  schemes  are  nearly  identical. 
Computational  details  of  the  non-centered  echmae 
ere  also  quits  different,  but  they  again  produce 
similar  steady  state  solutions. 

A  consistent  result  throughout  th* 
evaluations  was  that  the  non-centered  schemes 
produced  lower  stagnation  pressure  errors  than  did 
the  cantered  schemes.  Since  entropy  conservation  or 
stagnation  pressure  conservation  is  not  a 
conservation  property  of  th*  finite  difference 
equations,  wa  expect  thee*  errors  to  be  a  function 
of  the  truncation  errors  present  in  the  solution. 
The  truncation  error  appears  as  an  autrepy 
generation,  and  tha  stagnation  preasur*  error  is 
tha  entropy  error  raised  to  the  y/ly-D  power. 
Stagnation  pressure  loss  is  a  very  sensitive 
nee sure  of  the  entropy  generation  in  s  scheme. 

Tha  grid  si see  for  the  inviscid  problem*  were 
chosen  to  produce  smell  stagnation  pressure  errors 
for  tha  non-centered  schemes,  and  the  centered 
schemes  nearly  always  produced  larger  errors  on  the 
■asm  grid,  we  attribute  this  result  either  to  th* 
apaea  averaging  properties  of  the  solution 
operators  or  to  artificial  smoothing  to  eliminate 
double  sawtooth  errors. 

Accurate  Baler  flaw  simulations  are  possible 
with  the  non-centered  schemas  as  was  shown  by  th* 
supercritical  stator  axaapla.  Seth  methods 
reproduced  th*  shoek-fre*  solution  with  a  minimum 
of  effort,  and  both  produced  only  e  smell 
stagnation  pressure  loss  on  the  tsst  grid,  with  a 
locally  refined,  orthogonal  grid,  the  total  number 
of  grid  points  required  ean  be  reduced  to  a  wore 
manageable  mmber.  Our  ability  to  oouputa  accurate 
viscous  flow  solutions  is  increasing  rapidly.  Th* 
turbine  cases da  solution  is  fur  a  moderately 
difficult  geometry,  aad  tb*  computational  results 
are  reasonable.  Much  more  work  must  be  done  ia  this 
area  before  we  can  ooaputa  viscous  solutions  as 
accurately  as  for  inviscid  situations. 
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Figure  6a.  Surf  sen  Mach  number ,  Hi  bimp. 
Beam  Warming  method 
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Figure  6b.  Surface  Mach  nunbar,  Hi  bunp, 
modified  Beam-Warming  aachod 
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TASK  I  B:  ADVANCED  MEASUREMENTS  IN  TRANSONIC  COMPRESSORS 

Introduction 

This  research  effort  is  aimed  at  elucidating  design  point 
operation  loss  mechanisms  of  high  work,  high  mass  flow  transonic 
compressor  stages.  In  particular,  the  Air  Force  Aero  Propulsion 
Laboratory  High  Through  Flow  (HTF)  stage  is  being  examined  in  the 
M.I.T.  Blowdown  Compressor  Facility. 

The  M.I.T.  Gas  Turbine  Laboratory  has  developed  a  unique 
capability  for  measuring  the  time  dependent  three-dimensional  flow 
field  of  transonic  compressors.  Using  fast  response  semi-conductor 
diaphragm  pressure  transducers  mounted  on  traversing  probes,  the  time 
dependent  flow  fields  upstream  and  downstream  of  the  blade  rows  can 
be  mapped.  This  had  been  done  previously  for  the  MIT  stage  and  a 
NASA  20-inch  diameter  high  aspect  ratio  compressor  stage  and  is 
continuing  on  the  AFAPL  17-inch  diameter  HTF  stage. 

In  all  of  the  above  aerodynamic  measurements,  the  efficiency  of 
each  stage  was  obtained  by  direct  measurement  of  stagnation  pressure 
together  with  a  calculated  stagnation  temperature.  This  temperature 
was  found  using  the  Euler  equation  and  the  measured  tangential  Mach 
number.  This  calculation  rests  on  the  assumptions:  1)  that  the  flow 
is  steady  in  coordinates  rotating  with  the  rotor,  and  2)  that  there 
is  no  exchange  of  energy  between  stream-tubes  ( in  rotor  coordinates ) 
due  to  work  done  by  viscous  forces  or  due  to  heat  transfer.  Both  of 
these  assumptions  are  known  from  qualitative  observations  to  be  not 
completely  correct,  but  in  the  absence  of  a  direct  measurement  of  the 
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stagnation  temperature  at  the  rotor  exit  it  has  not  been  possible  to 
relax  the  assumptions  in  treating  the  details  of  the  rotor  outflow. 

For  this  reason,  we  have  started  work  on  a  technique  to  measure 
the  fluctuating  total  temperature  in  the  Blowdown  Compressor  as 
described  below. 

Time  Resolved  Total  Temperature  Probe;  Principle  of  Operation 

Time  resolved  measurement  of  gas  total  temperature  has  always 
been  a  difficult  task.  Conventional  thermocouples  lose  frequency 
response  above  1  KHZ.  Compensated  thermocouples  [1.1]  have  yet  to  be 
successfully  demonstrated.  Spectroscopic  techniques  using  flows 
seeded  with  NO  [1.2,  1.3],  or  2,3  butanedione  [1.4]  show  promise  as 
sensitive  non-instrusive  tools.  They  require,  however,  high  power, 
narrow  band  lasers  for  excitation  at  pulse  rates  sufficient  to  resolve 
the  flow  field.  For  transonic  compressors  where  significant  flow 
detail  occurs  at  frequencies  above  10  KHZ,  commercial  light  sources 
are  not  presently  available. 

We  have  decided  to  take  a  somewhat  simpler  approach  to  the 
problem  in  the  interests  of  returning  information  on  the  compressor 
flowfield  more  rapidly  than  might  result  from  the  more  elegant  but 
also  more  complicated  spectroscopic  techniques.  We  have  developed  a 
dual  wire  aspirating  probe  which  should  allow  measurement  of  both 
total  temperature  and  total  pressure  at  frequencies  above  10  KHZ. 
Conventional  hot  wires  measure  the  thermal  energy  lost  by  convection 
to  the  fluid.  They  are  thus  sensitive  to  the  wire/fluid  temperature 
difference  { Tw  -  rTip)  and  the  mass  flux,  pU,  by  the  wire,  where 
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Tw  is  the  wire  temperature,  TT,  the  fluid  total  temperature;  r 
the  recovery  factor;  p  the  fluid  density,  and  U  the  fluid  velocity. 
Since  the  voltage  from  a  constant  temperature  hot  wire  is 
proportional  to  the  square  root  of  the  power  dissipated  in  the  fluid. 

V*  =  f (pU) <TW  -  rTT)  ( 1 ) 

In  a  transonic  compressor,  all  the  flow  quantities  fluctuate 
temporally,  so  that  p,  U,  Tp,  and  Pip  (the  fluid  total  pressure)  are 
functions  of  time,  thus  making  conventional  hot  wire  data  exceedingly 
difficult  to  interpret  for  this  flowfield.  To  overcome  this  problem, 
we  have  placed  two  coplanar  hot  wires  in  the  1.5mm  diameter  channel 
of  an  aspirating  probe  (figure  1.1).  Since  the  convergent  exit  of  the 
channel  is  choked,  the  mass  flux  (pU)  is  now  a  function  only  of 
free  stream  total  temperature  (T-p)  and  total  pressure  (P-p), 
assuming  gas  composition  is  constant.  Thus,  equation  (1)  may  be 
written  as 


V1  =  f(PT//i^)  (Tw  -  rTip)  (2) 

In  order  to  verify  that  the  aspirating  probe  does  not  display 
the  usual  hot  wire  sensitivity  to  velocity  fluctuations,  the  response, 
of  the  probe  and  of  an  isolated  hot  wire,  to  the  vortex  shed  from 
a  right  circular  cylinder  were  compared.  When  the  cylinder  was 
unheated,  the  hot  wire  responded  to  the  velocity  fluctuations  but  the 
aspirating  probe  did  not  (figure  1.2).  When  the  wire  is  heated, 


-24- 


how  ever,  the  probe  responded  to  the  temperature  fluctuations,  thus 
confirming  the  performance  of  the  probe  to  first  order. 

Treating  the  hot  wire  as  a  right  circular  cylinder  and  adopting 
the  functional  form  of  the  heat  transfer  suggested  by  Collis  and 
Williams  [1.5].  Thus,  equation  2  becomes 

V*  =  [Ci  <PT//r^)n  +  C2]  (Tw  -  rT-j.)  (3) 

where  C^,  C2»  n  are  constants  to  be  determined  by  calibration 
(although  our  measured  values  are  in  good  agreement  with  those  of 
[1.5] ). 

Steady  State  Probe  Calibration 

As  for  any  instrumentation,  the  quality  of  the  data  can  be  no 
better  than  the  accuracy  of  the  calibration.  Initial  checks  of  the 
aspirating  probe  operation  were  performed  in  a  2.5cm  diameter  heated 
free  air  jet  in  which  the  total  pressure  (Mach  number)  and  the  total 
temperature  could  be  varied  independently.  The  range  of  Mach  number 
variation  is  0  to  1.0  and  of  total  temperature  20°C  to  150°C. 

Thus,  probe  calibration  in  air  could  be  established.  Also,  frequency 
response  was  roughly  estimated  by  using  the  vortex  shed  from  a  heated 
right  circular  cylinder  as  excitation.  This  excitation  frequency 
could  be  varied  up  to  about  8  KHZ  by  changing  the  velocity. 

The  initial  use  of  the  probe,  however,  is  in  the  M.I.T.  Blowdown 
Compressor  Tunnel  which  uses  an  argon- freon- 12  mixture  with  a  ratio  of 
specific  heats  of  1.4  as  the  working  fluid.  Thus,  a  test  rig  was 
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constructed  to  facilitate  calibration  with  this  gas  composition 
(figure  1.3).  The  test  rig  is  a  blowdown  facility  using  the  tanks  of 
the  blowdown  compressor  as  source  and  vacuum  dump.  The  supply  tank  is 
heated  by  an  atmospheric  pressure  steam  jacket  up  to  100°C.  The 
probe  is  tested  in  a  one  inch  diameter  channel  terminated  by  a  choked 
orifice.  The  Mach  number  in  the  channel  is  varied  by  changing  the 
orifice  diameter.  The  test  channel  is  heated  in  order  to  inhibit  the 
growth  of  a  thermal  bouundary  layer.  The  momentum  boundary  layer 
thickness  at  the  probe  position  is  approximately  10%  of  the  channel 
diameter.  The  data  which  follows  was  taken  at  a  channel  Mach  number 
of  0.19. 

The  calibration  procedure  is  to  take  data  by  varying  the  tank 
stagnation  pressure  at  constant  stagnation  temperature  while  operating 
at  a  particular  wire  overheat  ratio.  The  conditions  surveyed  include 
the  operating  range  of  the  Blowdown  Compressor  facility.  A  typical 
calibration  curve  is  presented  in  figure  1.4.  The  constants,  Ci, 

C-2,  and  n  are  determined  by  a  nonlinear  least  squares  data  fit. 

The  calibration  repeatability  appears  to  be  better  than  0.1%  in 
voltage  which  would  correspond  to  approximately  0.4%  in  adiabatic 
efficiency.  The  very  good  agreement  between  the  constants  derived 
from  calibration  runs  at  different  temperatures  implies  that 
interpolation  is  sufficient  for  any  intermediate  temperature. 


In  practice,  the  probe  data  is  reduced  with  an  iterative  routine 
using  the  experimental  calibration  curves.  This  routine  was  checked 
by  simultaneous  solution  of  equation  (3)  for  each  wire.  This  gave 
identical  results. 

As  a  further  check  on  probe  operation,  a  semiconductor  diaphragm 
(Kulite  XCQ-093)  total  pressure  probe  has  been  mounted  "piggyback"  on 
the  probe.  The  total  pressure  probe  has  the  same  inlet  diameter  as 
the  aspirating  probe  and  is  immediately  adjacent  to  it.  Thus,  total 
pressure  is  measured  in  two  ways.  To  date,  however,  the  fluctuating 
component  of  these  total  pressures  does  not  agree  well  (which  will  be 
discussed  later). 

Limits  to  Performance 

The  design  of  the  dual  wire  aspirating  probe  is  an  engineering 
compromise  between  spatial  resolution  (small  size)  and  hot  wire 
constraints  (wire  L/D);  and  angular  acceptance  and  frequency  response. 
The  hot  wire  length  to  diameter  ratio  must  be  kept  large  (above  200) 
in  order  for  thermal  end  effects  to  be  small.  Since  the  smallest 
practical  hot  wires  are  0.1  to  0.2  mils  in  diameter,  the  channel 
diameter  must  be  larger  than  1mm.  Because  of  anemometer  stability 
problems,  a  0.25  mil  diameter  wire  is  used  and  the  channel  diameter  is 
1.5mm.  This  is  also  consistant  with  construction  constraints.  In 
order  to  ensure  that  both  wires  exhibit  similar  sensitivity,  the  ends 
of  the  wires  are  copper  plated  so  as  to  equalize  the  active  length  of 


each  wire  to  within  1% 
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During  operation  in  the  blowdown  tunnel,  a  shift  in  D.C. 
calibration  of  0.2  «/•  0.5%  was  noted  after  each  run.  This  is  believed 
due  to  the  deposition  of  the  detonation  products  (mostly  fine  carbon) 
of  the  plastic  explosive  used  to  start  the  tunnel.  Careful  washing  of 
the  wires  in  distilled  water  after  each  run  returned  the  original 
calibration.  This  is  troublesome  since  most  wire  breakage  occurs 
during  handling.  A  future  solution  may  be  a  shutter  or  cap  on  the 
wall  cavity  in  which  the  probe  resides  during  the  explosion. 

The  basic  assumption  concerning  the  operation  of  the  probe  is 
that  both  wires  are  immersed  in  the  same  uniform  flowfield.  It  is 
only  to  the  degree  that  this  condition  holds  true  that  the  probe  can 
yield  meaningful  results.  This  assumption  cam  be  violated  in  two 
ways.  One  is  that  flowfield  contains  disturbances  whose  length  scales 
are  small  compared  to  the  probe  diameter.  This  is  certainly  true  but 
may  not  be  significant.  The  second  problem  would  arise  when  the  flow 
amgle  into  the  probe  is  large  enough  to  induce  separation.  One  way  of 
improving  the  probe's  response  to  off  amgle  amd  small  scale  dis¬ 
turbances  would  be  to  increase  the  L/D  of  the  channel  before  the 
wires,  thus  mixing  out  the  non -uniformities.  This  approach  cannot  be 
taken,  however,  because  the  growth  of  the  resultamt  thermal  boundary 
layer  would  damp  out  the  high  frequency  components  of  interest.  The 
present  design  was  chosen  to  keep  the  thermal  boundary  layer  height  at 
the  wire  plame  below  10%  of  the  chamnel  diameter.  The  off-amgle 
response  of  the  probe  has  been  checked  to  ±20°  with  no  deviation 


noted 


Another  assumption  in  the  data  analysis  is  that  the  channel  flow 
in  the  probe  behaves  in  a  steady  state  quasi- ID  manner.  This  was 
checked  by  examining  the  solution  to  the  time  dependent  continuity  and 
momentum  equations  for  the  conditions  relevant  to  the  probe's  opera¬ 
tion.  No  significant  deviation  from  the  steady  state  assumption  was 
apparent. 

An  additional  limit  to  performance  is  the  basic  signal  to  noise 
ratio  of  each  hot  wire.  At  the  present  low  total  pressure  required 
for  steady  state  (i.e.  constant  corrected  speed  and  weight  flow) 
operation  in  the  blowdown  tunnel,  the  blade  wake  temperature  excess  is 
only  four  times  the  noise  level.  For  this  reason,  the  temperature 
derived  from  using  the  piggyback  Kulite  total  pressure  probe  (which 
has  relatively  little  noise)  and  any  one  of  the  two  wires  gives  a  more 
accurate  reading  than  the  temperature  derived  from  both  wires 
together.  This  problem  was  recognized  after  the  data  presented  herein 
was  taken  and  will  be  improved  in  a  later  data  set.  Also,  the  inlet 
total  pressure  will  be  doubled  with  the  flywheel  installed,  thus  also 
increasing  the  signal  to  noise  of  the  probe. 

The  final  constraint  on  confirmation  of  overall  probe  accuracy  is 
the  lack  of  an  accurate  high  frequency  calibration  of  the  aspirating 
probe.  The  shed  vortex  method  is  qualitative  at  best  since  Reynolds 
number  constraints  limit  the  size  of  cylinder  and  thus  the  wake. 
Although  the  time  dependent  continuity  analysis  indicates  that  channel 
flow  acts  as  though  in  steady  state,  experimental  verification  is 
strongly  desirable.  Unfortunately,  it  is  extremely  difficult  to 
generate  a  well-characterized  strongly  fluctuating  flow  field  at 


frequencies  upward  of  1  KHZ.  We  are  presently  considering  several 


methods  of  accomplishing  this  including  a  two-foot  diameter  chopper 
disk  and  a  high  speed  "lawn  sprinkler"  concept.  In  both  cases,  these 
are  complex,  high  speed  rotating  devices  requiring  careful  mechanical 
and  aerothermal  design.  Plans  are  now  to  construct  such  a  device 
during  the  next  year. 

Preliminary  Measurements 

The  aspirating  probe  has  been  used  to  survey  the  outflow  of  the 
AFAPL  High  Through  Flow  (HTF)  rotor.  The  rotor  is  operated  with  the 
probe  inserted  in  a  stator  gap  directly  behind  the  rotor.  The  probes 
are  traversed  from  tip  to  hub  in  15  to  20  ms.  These  measurements  were 
all  made  before  the  installation  of  the  flywheel  on  the  HTF  rotor. 
Thus,  the  corrected  speed  changes  approximately  1.5%  during  the  test 
time. 

The  spanwise  distribution  of  total  pressure  ratio  across  the 
rotor  is  presented  in  figure  1.5.  The  total  pressure  probe  is  the 
semiconductor  unit  mounted  piggyback  on  the  aspirating  probe.  The 
four-way  probe  is  a  high  frequency  cylinder  probe  illustrated  in 
figure  1.6.  The  AFAPL  stator  L.E.  refers  to  conventional  steady  state 
pitot  probes  located  along  the  stator  leading  edge.  These  are  conven¬ 
tional  measurements  taken  at  AFAPL.  For  the  two  high  frequency 
measurements,  the  data  are  time  averages.  Note  that  the  agreement 
among  the  measurements  is  not  very  close.  This  difference  is  of  the 
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order  of  3-5%  at  a  given  spanwise  location.  Discrepancies  between 
steady  state  and  time  resolved  data  have  been  noted  before  [1.6] 
(although  not  necessarily  explained). 

Evidence  from  ongoing  testing  at  NASA  Lewis  Research  Center 
suggests  that  the  spanwise  pressure  distribution  of  a  transonic  com¬ 
pressor  may  continuously  vary  to  some  degree,  even  though  the  average 
operating  point  is  steady.  To  explore  this,  surveys  using  the  total 
pressure  probe  during  three  nominally  identical  runs  is  shown  in 
figure  1.7.  Although  the  curves  have  generally  the  same  shape  and 
average  to  the  same  total  pressure,  the  particular  value  at  any 
spanwise  location  may  vary  by  as  much  as  3%  from  test  to  test. 

Another  aspect  of  this  variation  can  be  seen  in  figure  1.8  which 
shows  the  wall  static  pressure  approximately  0.2  chord  upstream  of  the 
tip.  Here,  the  waveforms  have  the  same  overall  shape  but  are  clearly 
not  identical  at  the  few  percent  level.  (The  traces  are  displaced 
vertically  for  clarity. )  This  is  somewhat  surprising  since  the  blow¬ 
down  tunnel  should  be  an  extremely  repeatable  facility  for  compressor 
testing  (even  the  small  acoustic  disturbances  in  the  supply  tank  are 
identical  from  test  to  test) .  Even  if  time  resolved  transonic 
compressor  behavior  is  unsteady,  it  should  be  repeatable  in  the  blow¬ 
down  facility.  It  is  suggested  that  the  variation  observed  in  these 
three  tests  is  due  to  differences  in  the  initial  angular  position  of 
the  rotor  at  the  instant  the  test  starts.  Thus,  the  flow  variations 
would  be  a  function  of  three  variables  —  blade-to-blade  variations 
(and  thus  rotor  angular  position),  inflow  perturbations,  and  any 
inherent  rotor  unsteadyness.  Exactly  in  what  manner  these  three 


i 


phenomena  interact  is  not  well  understood  at  this  time  but  is  a  matter 
of  continuing  research.  In  an  attempt  to  understand  the  causes  of  the 
test-to-test  variation  in  the  blowdown  compressor,  future  tests  will 
be  started  with  the  rotor  at  identical  (to  1/40  blade  spacing)  angular 
positions.  It  is  believed  that  this  should  substantially  reduce  this 
variation  (although  not  the  unsteadiness). 

With  the  level  of  test-to-test  variation  in  mind,  we  can  now 
compare  the  measurements  from  different  probes  at  the  same  location 
during  different  tests.  Time  resolved  measurements  of  rotor  outflow 
total  pressure  is  presented  in  figure  1.9  for  probes  of  differing 
configurations  and  size.  The  sphere  probe  is  approximately  5mm  in 
diameter;  the  four-way  cylinder  probe,  3mm;  the  total  pressure  probe 
with  a  Kulite  XCQ-093  tranducer,  3mm;  and  the  total  pressure  probe 
with  a  Kulite  XCQ-063,  2mm.  The  D.C.  levels  are  the  same  but  the 
traces  have  been  displaced  vertically  for  clarity.  As  would  be 
expected,  the  frequency  response  improves  as  the  probe  size  is 
reduced.  More  importantly,  however,  the  symmetry  of  the  waveform 
changes,  especially  near  the  tip,  as  the  probe  size  is  reduced.  Thus, 
conclusions  as  to  blade -to- blade  flow  details  must  be  made  with 
extreme  caution.  The  time  average  of  these  signals  vary  no  more  than 
the  test-to-test  differences  shown  in  figure  1.7.  Further  exploration 
of  the  probe  size  effect  will  include  construction  of  a  1.2mm  diameter 
high  frequency  total  pressure  probe.  Also,  any  high  frequency  cali¬ 
bration  facility  constructed  for  use  with  the  aspirating  probe  will 
also  be  utilized  to  check  the  other  aerodynamic  probes. 
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Comparison  of  the  spanwise  distribution  of  the  time  average  of 
the  total  temperature  a)  measured  with  the  aspirating  probe  and  b) 
inferred  using  the  Euler  Turbine  Eq.  from  angle  measurements  made  with 
the  four-way  cylinder  probe  is  shown  in  figure  1.10.  The  resultant 
adiabatic  efficiency  is  shown  in  figure  1.11.  It  is  important  to  note 
that  these  averages  differ  no  more  than  the  test-to-test  flow  varia¬ 
tion.  Thus,  on  a  time  averaged  basis,  one  may  conclude  that  the  Euler 
Equation  is  an  accurate  approximation  of  the  rotor  flow. 

This  is  not  necessarily  the  case  for  the  blade-to-blade  flow 
detail  shown  in  figures  1.12-1.14.  Here,  ten  blade  ensemble  averages 
of  the  flew  near  the  nub,  midspan,  and  tip  are  presented  along  with 
their  respective  power  spectral  densities.  In  this  case,  the 
difference  in  the  blade-to-blade  variation  is  significant  even  though 
the  average  flow  quantities  are  equal.  In  the  case  of  the  tip 
measurement  (R/Rt  =  0.95)  much  of  the  variation  is  due  to 
differences  in  the  total  pressure  waveform,  perhaps  due  to  probe  size 
effects  mentioned  earlier.  This  is  not  the  case  for  the  inner 
stations,  however,  where  the  variation  i‘!s  principally  in  the  total 
temperatures.  Thus,  the  Euler  Equation  may  not  be  an  accurate  repre¬ 
sentation  of  the  time  resolved  flow  field.  Further  analysis  work  is 


planned  in  this  area 
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Plans  for  Future  Work 

Efforts  for  the  next  reporting  period  include  both  further 
testing  and  analysis.  The  HTF  rotor  will  be  retested  with  the 
recently  acquired  flywheel  installed.  This  will  reduce  the  corrected 
speed  and  weight  flow  variations  to  less  than  0.5%  and  also  increase 
the  Reynolds  number  by  a  factor  of  two. 

A  high  frequency  calibration  facility  will  be  designed  and 
constructed  to  test  the  frequency  response  of  the  aspirating  and  other 
probes.  Also,  work  will  proceed  on  improving  the  signal-to-noise 
ratio  of  the  aspirating  probe.  Analysis  on  the  data  will  continue 
with  emphasis  placed  on  understanding  the  implications  of  the  time 
resolved  temperature  data. 
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Fig.  1.2a  —  Probe  vs.  hot  wire  response  in  a  vortex  sheet 
shed  by  an  unheated  cylinder. 
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Scale  Drawing  of  Aspirating  Probe  Calibration  Facility 
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Fig.  1.5  —  Comparison  of  time  averages  of  4-way  cylinder 

probe  and  total  pressure  vs.  pitot  measurements 
at  AFAPL. 
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4-way  high  frequency  cylinder  probe 
(che  light  spots  are  transducers). 
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Fig.  1.7  —  Variation  in  radial  distribution  of  stagnation  pressure 


a*  aeasured  by  total  pressure  probe  during  three 
nominally  identical  tests. 
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Comparison  of  wall  static  pressure  measured  0.2  chords  upstream  of 
rotor  leading  edge  during  3  nominally  identical  runs.  The  curves 
have  been  displaced  vertically  for  clarity. 
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Variation  in  measured  total  pressure  for  4  different  probes  during 
a  tip  to  hub  traverse  during  4  different  tests.  The  traces  are 
displaced  vertically  for  clarity. 
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Fig.  1.10  —  A  Comparison  of  time  averaged  total  temperature 
measured  with  the  aspirating  probe,  conventional 
thermographs  at  AFAPL,  and  deduced  from  angle 
measurements  using  Euler's  Equation. 
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1.12a  —  Cosnparsion  of  total  pressure  ratio,  total 
temperature  ratio  and  efficiency  between 
4  Way  Probe  and  Aspirating  Probe  at  R/Rt».95 
10  cycle  ensemble  averaged. 
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Fig.  1.12b  — 


Frequency  spectrum  of  total  pressure  ratio,  total 
temperature  ratio  and  adiabatic  efficiency  from 
4  way  probe  &  aspirating  probe  at  R/Rt«.95; 

10  cycle  ensemble  averaged. 
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Fig.  1.13a  —  Comparison  of  total  pressure  ratio,  total 
temperature  ratio  and  efficiency  between 
4  Way  Probe  and  Aspirating  Probe  at  R/Rt- . 88 
10  cycle  ensemble  averaged. 
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Fig.  1.13b  —  Frequency  spectrum  of  total  pressure  ratio,  total 
temperature  ratio  and  adiabatic  efficiency  from 
4  way  probe  &  aspirating  probe  ar  R/Rt*.88; 

10  cycle  ensemble  averaged. 
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Fig.  1.14b  — 


Frequency  spectrum  of  total  pressure  ratio,  total 
temperature  ratio  and  adiabatic  efficiency  from 
4  way  probe  &  aspirating  probe  ar  R/Rt».73; 

10  cycle  ensemble  averaged. 
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task  II:  MECHANISMS  OF  COMPRESSOR  STABILTY  ENHANCEMENT 
USING  CASING/HUB  TREATMENT 

During  the  past  year  we  have  tested  the  hub  treatment  using  a 
high  stagger  stator  configuration,  in  order  to  promote  a  wall  stall. 
The  result  was  a  substantial  improvement  in  stator  stall  margin  and 
performance.  This  is  the  first  time  that  such  improvements  have  been 
documented . 

We  have  also  continued  efforts  on  the  modelling  of  the  inter¬ 
action  between  the  flow  in  the  grooves  and  that  in  the  endwall 
(although  the  main  emphasis  has  been  on  the  experiments).  The 
detailed  results  are  given  in  the  M.S.  Thesis  of  Mr.  P.  Cheng. 

However,  a  paper  has  also  been  written  which  summarizes  this  work,  and 


this  is  presented  below 


EFFECTS  OF  COMPRESSOR  HUB  TREATMENT 
ON  STATOR  STALL  AND  PRESSURE  RISE 


P.  Cheng,  M.E.  Prell,  E.M.  Greitzer,  C.S.  Tan 
Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


ABSTRACT 


An  investigation  has  been  conducted  of  the  changes,  in  stator 
stall  margin  and  performance,  which  occur  due  to  a  slotted  hub 
treatment  rotating  beneath  an  axial  compressor  stator  row.  Two 
different  stator  rows,  one  with  low  stagger  and  one  with  high  stagger, 
were  tested  both  with  solid  wall  and  with  hub  treatment.  The 
experiments  showed,  for  the  first  time,  that  (with  the  high  stagger 
blading)  the  hub  treatment  was  very  effective  in  delaying  the  outset 
of  stator  stall  and  in  increasing  the  peak  stator  static  pressure 
rise.  For  this  blading,  measurements  of  the  stator  exit  flow  field 
appear  to  indicate  that  the  inception  of  stall  was  associated  with  a 
wall  stall.  With  the  treatment  the  large  blockage  in  the  hub  region 
was  greatly  decreased  and,  near  the  hub,  the  total  pressure  at  the 
stator  exit  was  found  to  be  higher  than  that  at  stator  inlet.  These 
results,  in  accompaniment  with  the  data  taken  using  the  low  stagger 
stator,  support  the  hypothesis  that  endwall  treatment  is  effective 
when  the  type  of  stall  that  occurs  is  wall  stall. 
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Introduction 

It  has  been  known  for  some  years  that  the  application  of  grooves  or 
slots  over  the  rotor  tips  in  the  casing  of  an  axial  flow  compressor  can 
have  a  large  effect  on  compressor  stall  margin  and  performance.  Experi¬ 
mental  observations  (1,  2)  have  also  shown  that  casing  treatment  can 
substantially  reduce  the  endwall  boundary  layer  blockage  in  the  rotor  tip 
region.  It  has  been  hypothesized  that  this  blockage  reduction  results 
from  strong  flow  injection  out  of  the  treatment,  and  that  the  low-momen¬ 
tum  flow  in  the  endwall  region  is  energized  due  to  this  injection  of  high 
velocity  fluid. 

It  has  also  been  found  that  one  can  often  make  a  simple  but  useful 
distinction  between  two  types  of  compressor  stall.  One  of  these,  termed 
"blade  stall"  is  a  roughly  two-dimensional  type  of  stall  where  a 
significant  portion  of  the  blade  span  has  a  large  wake  due  to  substantial 
thickening  or  separation  of  the  suction  surface  boundary  layer.  The 
other,  termed  "wall  stall,"  is  associated  with  the  endwall  boundary  layer 
separation  (1,  2),  and  it  appears  that  casing  treatment  will  only  be 
effective  when  the  type  of  stall  is  wall  stall  (1,  2). 

A  third  feature  which  has  emerged  not  only  from  (1,  2)  but  from  flow 
visualization  reported  in  (3),  is  that  relative  motion  between  endwall 
and  blade  is  an  important  element  for  the  success  of  the  grooves  in 
reducing  the  blockage.  In  view  of  this,  it  is  natural  to  ask  whether  a 
rotating  "hub  treatment"  below  a  row  of  cantilevered  stator  blades  would 
also  be  effective  in  improving  the  stall  margin.  If  this  is  so,  the 
geoaetry  would  also  be  more  amenable  to  investigation  of  the  details  of 
flow  in  the  blade  passage,  as  well  as  provide  a  different  perspective 


from  which  to  view  the  fluid  dynamic  phenomena 
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In  the  present  work,  the  effect  of  hub  treatment  on  the  stall 
margin  and  performance  of  a  "hub  critical”  stator  was  examined.  Two 
sets  of  blading  configurations  were  tested,  each  with  a  solid  wall 
and  with  hub  treatment.  It  was  found  that  for  the  first  set  of 
blading  (low  stagger)  the  application  of  hub  treatment  did  not  have  a 
major  effect  on  the  stator  performance.  However,  for  the  second  set 
of  blading  (high  stagger)  the  hub  treatment  improved  the  compressor 
performance,  delaying  the  inception  of  stall  and  markedly  increasing 
the  peak  stator  static  pressure  rise. 

Experiment  Design 

It  is  fundamental  in  the  design  of  this  experiment  to  ensure  that 
the  stator  hub  endwall  is  stall  limiting.  The  design  used  is  aimed  at 
having  the  rotor  (and  IGV)  act  only  as  a  "flow  generator"  for  the 
stator.  Thus,  in  choosing  the  blade  setting  angles,  it  is  desired 
that  the  rotor  has  a  low  loading  relative  to  the  stator.  The 
requirement  of  high  stator  hub  loading  can  be  satisfied  by  using  a 
rotor  blade  with  low  twist.  This  creates  a  low  total  pressure  rise 
across  the  rotor  hub  relative  to  the  tip,  which  loads  the  stator  hub. 

When  the  initial  experiments  were  carried  out  in  the  preesent 
facility,  a  severe  restriction  was  that  the  stator  discharge  was  to 
ambient  pressure.  Hence,  as  described  in  (4),  we  were  forced  to  use 
blading  of  rather  low  stagger,  which  tends  not  to  be  conducive  to  wall 
stall.  In  accordance  with  this  (as  described  in  detail  in  (4))  the 
type  of  stall  encountered  in  these  initial  experiments  appeared  to  be 
blade  stall,  rather  than  pure  wall  stall. 

During  the  second  set  of  experiments,  which  were  carried  out  with 
an  exhaust  fan,  however,  we  had  the  freedom  to  increase  the  blade 
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stagger;  this  tended  to  promote  a  wall  stall.  In  addition ,  in  order 
to  achieve  a  wall  stall,  based  on  the  correlation  presented  in  (5), 
the  stator  hub  clearance  was  doubled  to  four  percent  of  the  blade 
chord  (or  4.5  percent  of  the  staggered  gap  based  on  exit  angle.)  A 
summary  of  the  two  stator  configurations  used  is  presented  in  Table  1. 
In  order  to  distinguish  these  two,  we  will  refer  to  them  below  as  the 
low  stagger  and  high  stagger  stators  respectively. 


TABLE  1 

STATOR 

BLADING  GEOMETRY 

Low  Stagger 

High  Stagger 

Midspan  stagger  angle 
(measured  from  axial) 

22.5* 

42. 5» 

Camber 

30. 0» 

30.0° 

Solidity  at  midspan 

1.0 

1.0 

Blade  clearance 

0.8mm 

1.5mm 

Number  of  blades 

45 

45 

Chord 

38.0mm 

38.0mm 

EXPERIMENTAL  FACILITY 

The  experiment  was  conducted  on  a  single  stage  research 
compressor  driven  by  a  variable  speed  D.C.  motor.  This  is  a  0.59m 
diameter  single  stage  compressor  with  a  hub/tip  ratio  of  0.75.  It  was 
modified  for  casing  treatment  as  well  as  hub  treatment,  and  the  former 
was  used  throughout  all  the  tests,  as  this  tended  to  increase  the  hub 
loading.  A  cross  sectional  schematic  of  the  compressor  showing  the 
blading  and  treatment  locations  is  given  in  Figure  1. 

Most  of  the  measurements  were  taken  at  a  blade  speed  of  75m/s  (at 
the  mean  radius)  for  the  low  stagger  stator  and  at  70m/s  for  the  high 
stagger  stator.  Runs  were  also  made  at  speed  from  55m/s  to  95m/s  to 
check  Reynolds  number  dependence  (which  was  found  to  be  of  only  slight 
influence  on  the  observed  changes  due  to  the  hub  treatment ) •  The 
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Reynolds  number  based  on  blade  chord  at  the  stator  midspan  (at  the 
stall  point)  was  nominally  1.5  x  10s  for  the  first  set  of  tests  and 
1.0  x  10s  for  the  second. 

The  pressure  instrumentation  consisted  of  20  total  pressure 
"Kiel"  probes  and  24  hub  and  casing  static  pressure  taps.  Acquisition 
of  the  pressure  data  was  done  by  a  48  channel  scanivalve  pressure 
scanner  operated  by  a  digital  microcomputer.  A  self-nulling  radial 
traverse  mechanism  was  used  to  determine  radial  profiles  of  total 
pressure  and  flow  angle  and  a  hot  wire  anemometer  was  used  to  enable 
accurate  determination  of  the  stall  point. 

The  data  to  be  presented  is  non-dimensionalized  in  a  standard 
fashion:  velocities  are  non-dimensionalized  by  U,  the  mean  blade 
speed,  and  pressure  difference  are  non-dimensionalized  by  either 
1/2 pU*  or  by  the  relevant  stator  inlet  dynamic  pressure.  Other 
quantities  used  for  non-dimensionalization  or  for  reference  levels 
will  be  defined  as  they  are  introduced.  In  all  the  figures  "SW"  is 
used  to  denote  the  solid  wall  and  "HT"  is  used  for  hub  treatment. 
Further  details  of  the  facility  are  given  in  (4)  and  (6). 

TREATMENT  DESIGH 

The  hub  treatment  consists  of  axial  slots  skewed  at  a  60°  angle 
to  the  radial  direction.  These  slots  have  been  found  to  provide  large 
increases  in  stall  margin  when  used  as  a  rotor  casing  treatment.  The 
slots  rotate  under  the  middle  90%  of  the  stator  for  the  high  stagger 
rotor  and  under  the  middle  70%  for  the  low  stagger.  The  slot  spacing 
is  such  that  the  slotted  area  is  twice  that  of  the  solid  (land)  area. 
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The  slot  aspect  ratio  (axial  length/tangential  width)  is  2.0,  and  the 
radial  depth  is  30%  of  the  axial  length.  The  actual  dimensions  and 
treatment  geometry  are  shown  in  Figure  2.  The  hub  treatment  was 
constructed  by  machining  a  wide  circumferential  groove  in  the  cylinder 
and  glueing  pre-machined  plexiglass  plates  into  the  groove. 

EXPERIMENTAL  RESULTS 

In  this  section  we  first  show  the  effects  of  the  hub  treatment  on 
stall  margin  and  compressor  speedline  characteristics,  then  compare 
the  untreated  and  treated  traverse  data  and  discuss  these  results. 

Low  stagger  stator 

For  the  low  stagger  blading,  the  application  of  hub  treatment  to 
the  compressor  described  above  resulted  in  only  a  small  reduction  of 
the  stalling  flow  coefficient1  (approximately  2%)  and  an  increase  in 
the  stator  pressure  rise  characteristic  of  approximately  5%  at  the 
peak.  The  stator  pressure  rise  characteristic  from  this  test,  as 
determined  from  the  casing  wall  static  taps  just  ahead  of  and  just 
behind  the  stator,  is  shown  in  Figure  3.  The  vertical  axis  is  nondimen- 
sionalized  wall  static  pressure  rise  across  the  stator  and  the  hori¬ 
zontal  axis  is  the  axial  velocity  paramter  (C^/V)  based  on  mean  inlet 
axial  velocity.  The  stall  points,  i.e.,  the  points  of  inception  of 
rotation  stall  in  the  stator  are  also  indicated.  (For  reference  the 
rotor  stall  point  was  far  below  these  at  Cc/U  *  0.35.)  Rotating  stall 

*The  point  at  which  rotating  stall  was  encountered. 


onset  was  clearly  detectable  on  an  oscilloscope  vhich  displayed 
the  response  of  hot  wire  survey  at  stator  exit  20%  span  away  from  hub 


through  an  anemometer. 

Although  there  is  an  increased  static  pressure  rise  across  the 
stator,  due  to  the  hub  treatment,  at  all  flow  coefficients  tested,  the 
increase  is  much  less  than  the  large  increases  that  have  been 
experienced  with  casing  treatment.  One  additional  aspect  described  in 
(4)  is  that  there  was  a  decrease  of  stator  deviation  angle  and 
blockage  with  treatment  compared  to  the  untreated  case,  and  presumably 
this  improvement  is  the  reason  for  the  increase  in  treated  3tator 
static  pressure  rise.  However,  from  the  analysis  of  stator  wakes,  it 
was  seen  that  a  thick  wake  was  present  from  hub  to  midspan  both  with 
treatment  and  without.  This  data,  as  well  as  other  supporting 
evidence  given  in  (4),  implied  that  the  inception  of  stall  was 
associated  with  a  so-called  "blade  stall"  type  of  separation. 

High  stagger  stator 

For  the  high  stagger  blading  the  application  of  hub  treatment  of 
the  compressor  resulted  in  a  much  more  substantial  change.  The  flow 
coefficient  at  stall  onset  was  decreased  by  12%  and,  more 
impressively,  the  peak  stator  static  pressure  rise  was  increased  by 
over  50%.  Figure  4  shows  a  comparison  of  the  overall  total  to  static 
pressure  rise  characteristic  [(pex  -  ptin)/(pU  /2)  vs.  (C^/U)] 
for  the  complete  stage,  with  and  without  the  hub  treatment.  The 
stator  stall  points,  also  indicated  here,  are  at  Cx/U  **  .315  and  .275 
(for  the  solid  wall  and  hub  treatment  respectively),  with  an  uncer- 


tainty  of  approximately  one  percent.  Note  that  although  the  C*/U 
values  are  somewhat  lower  than  normal,  it  should  be  emphasized  that 
this  is  due  to  the  high  rotor  stagger,  and  the  stagger  and  flow  angles 
associated  with  the  stator  are  representative  of  those  encountered  in 
modem  compressors.  Due  to  the  large  pressure  rise  in  the  rotor,  the 
improvement  in  stator  performance  is  not  really  brought  out  in  this 
figure.  However,  the  figure  does  show  that  the  overall  stage 
characteristic  is  negative  and  that  the  rotor  is  operating  far  from 
stall,  as  should  be  the  case  in  this  experiment. 

A  comparison  of  the  stator  static  pressure  rise  characteristic  is 
shown  in  Figure  5,  which  shows  APstator/( PU*/2 )  plotted  versus 
C^/U  for  the  solid  wall  (SW)  and  the  hub  treatment  (HT).  It  is 
evident  that  the  inception  of  stall  has  been  delayed  and  the  static 
pressure  rise  has  been  markedly  increased  due  to  the  hub  treatment. 

These  stator  speed  lines  were  checked  several  times  with  and 
without  treatment.  It  was  found  that  there  was  some  slight 
deterioration  in  performance  over  long  test  times  (due  presumably  to 
the  accumulation  of  dirt,  which  could  be  seen  to  build  up  on  the 
compressor  blades).  However,  the  maximum  difference  found  due  to 
this  effect  was  a  flow  coefficient  shift,  with  both  solid  wall  and 
hub  treatment,  of  2%  or  less,  i.e.,  the  changes  were  an  order  of 
magnitude  less  than  the  differences  between  the  solid  wall  and  hub 


treatment  configurations 
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Radial  Traverse  Total  Pressure  Data 

In  order  to  examine  the  feature  of  the  flow  in  more  detail, 
radial  traverses  were  taken  with  the  solid  wall  and  with  hub  treatment 
at  a  "solid  wall  near  stall  point",  C^/V  =  0.33.  A  comparison  of 
the  (arithmetic)  pitch  average  total  pressure  distribution  at  the 
rotor  and  stator  exit  for  hub  treatment  and  solid  wall  at  this  near 
stall  point  is  presented  in  Figure  6.  The  symbols  RTE  and  STE  refer 
to  the  locations  downstream  of  the  rotor  and  stator  respectively,  as 
shown  in  Figure  1,  and  again,  SW  denotes  solid  wall  and  HT  denotes  hub 
treatment.  In  the  figure  the  vertical  axis  is  radial  location 
(percent  span)  and  the  horizontal  axis  is  the  non-dimensionalized 
total  pressure.  The  reference  pressure  is  ambient  (compressor  inlet 
total  pressure)  and  the  dynamic  pressure  used  for 

non-dimensionalization  is  based  on  the  midspan  stator  inlet  velocity. 
The  reason  for  this  latter  convention  is  that  due  to  the  very 
non-uniform  total  pressure  profile,  there  is  no  one  "freestream"  value 
of  the  stator  inlet  dynamic  pressure;  however,  the  midspan  quantity  is 
a  convenient  representative  value. 

In  examining  the  stator  exit  traverses  it  should  be  mentioned 
that  there  was  an  axial  gap  of  3mm  behind  the  stator  (see  Figure  1) 
between  the  rotating  hub  piece  and  the  stationary  annulus,  with  a 
dead-ended  cavity  below.  Although  this  may  affect  the  absolute  values 
of  total  pressure  and  flow  in  the  hub,  it  is  to  be  emphasized  that  it 
is  the  comparisons  between  solid  wall  and  hub  treatment  that  are  of 


most  interest  here,  and  the  gap  geometry  is  the  same  for  both  these 
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cases.  The  comparisons  should  thus  be  valid,  especially  since  the 
differences  are  quite  substantial. 

It  can  be  seen  that  the  rotor  exit  total  pressure  profiles  are 
nearly  identical  so  that  for  both  smooth  wall  (SW)  and  hub  treatment 
(HT)  the  (stator)  inlet  conditions  are  the  same.  However,  the  two 
stator  exit  profiles  are  quite  different.  With  the  solid  wall  there 
is  a  a  large  total  pressure  defect  (and  also  a  large  blockage)  near 
the  stator  hub.  In  contrast,  for  the  treated  case  the  stator  exit 
total  pressure  is  higher  than  the  rotor  exit  total  pressure  near  the 
hub.  It  is  evident  that  work  has  been  done  on  the  flow  in  the  endwall 
region  due  to  the  fluid  injected  from  the  grooves  and  that  the  heavy 
losses  (and  blockages)  in  the  hub  region,  have  been  reduced.  (Note 
that  for  this  set  of  experiments  the  flow  is  essentially  incom¬ 
pressible  so  there  Is  no  net  flow  out  of  the  treatment.  However, 
due  to  the  fact  that  the  flow  entering  and  leaving  the  grooves  has  a 
very  different  velocity,  there  can  be  a  net  momentum  flux  out  of  the 
grooves. ) 

Another  feature  seen  in  Figure  6  is  that  even  at  midspan  the 
treated  stator  exit  total  pressure  with  hub  treatment  is  apparently 
higher  than  that  which  occurs  with  the  solid  wall.  However,  one  would 
not  expect  the  mixing  effects  (i.e.  shear  forces)  to  extend  to  this 
radial  station  from  the  endwall  region.  In  fact  the  reason  for  this 
total  pressure  change  is  that  there  is  a  large  blockage  at  the  hub 
region  with  the  solid  wall,  and  the  flow  near  the  hub  is  forced  to 
move  outwards.  In  consequence,  there  is  a  considerable  radial 
streamline  shift  that  occurs  between  the  RTE  to  STB  axial  locations 
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for  this  configuration.  This  radial  shift  (convection)  of  the 
non-uniform  total  pressure  profile  causes  the  apparent  total  pressure 
change  in  the  midspan  region,  as  will  be  discussed  below. 

Flow  Angle 

Figure  7  presents  comparison  of  the  absolute  flow  angle  a, 
measured  from  the  axial  direction.  The  figure  shows  the  midpitch 
radial  distribution  of  angles  with  smooth  wall  and  with  hub  treatment 
at  Cx/U  =0.33.  It  is  seen  that  again  the  two  rotor  exit  profiles 
are  very  similar,  the  difference  between  the  two  being  within  0.5°  at 
all  radii.  In  contrast,  the  two  stator  exit  profiles  are  quite 
different.  With  the  treatment  there  is  an  increase  in  angle  very  near 
the  hub  but  a  decrease  in  angle,  i.e.,  an  increase  in  turning,  in  the 
midspan  region  due  to  the  treatment?  similar  results  were  obtained  in 
the  low  stagger  case. 


Circumferential  Traverses 

Another  way  to  examine  the  effect  of  treatment  on  the  flow  field 
is  to  compare  the  circumferential  distribution  of  total  pressure  at 
the  STE  station  for  the  solid  wall  and  the  hub  treatment.  This 
comparison  is  made  in  Figure  8,  which  shows  plots  of  total  pressure 
difference  across  the  stator  blade  row  divided  by  local  inlet  dynamic 
head,  versus  circumferential  position,  for  several  different  span 
locations.  The  total  pressure  used  for  reference  in  this  figure  is 
the  rotor  exit  total  pressure,  which  has  been  converted  to  the 
appropriate  stator  exit  radial  position  using  a  streamline  curvature 
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data  reduction  program  (7).  The  dynamic  pressure  used  for 
nondimensionalization  is  based  upon  the  local  stator  inlet  velocity. 

Although  the  reference  pressures  at  inlet  are  slightly 
different  due  to  the  (radial)  streamline  shifts,  it  can  be  seen  that 
there  are  major  differences  between  these  two  cases  at  the  region  6% 
span  away  from  hub.  For  the  solid  wall  the  wake  is  very  deep  with 
large  loss.  With  hub  treatment  the  wake  profile  is  flatter  and,  more 
importantly,  the  total  pressure  is  higher  than  the  rotor  inlet  total 
pressure  at  all  circumferential  positions. 

Further  evidence  of  the  differences  between  solid  wall  and  hub 
treatment  is  seen  in  the  comparison  of  the  wake  distribution  at  the 
26%  span  location.  With  the  solid  wall  the  wakes  are  fairly  thin; 
this  is  viewed  as  another  indication  that  the  onset  of  stall  is  due  to 
a  pure  wall  stall  at  the  hub.  With  the  hub  treatment,  however,  there 
is  a  "spike"  in  total  pressure  between  near  the  pressure  surface  side. 

We  take  this  to  imply  that  this  is  an  indicatioin  of  a  jet  of  high 

total  pressure  fluid.  This  is  consistent  with  a  picture  of  strong  ; 

flow  injection  out  of  the  slots  impingent  on  the  pressure  surface,  and  I 

I 

deflection  rearwards.  At  the  midspan  region  and  outboard  the  wake  ; 

distributions  are  similar  for  the  solid  wall  and  the  treatment, 

although  the  magnitudes  are  slightly  different. 

The  effect  of  hub  treatment  on  the  total  pressure  "loss"  profile  , 

(based  on  the  axisymmetric  data  analysis  program)  is  shown  in  Figure 

I 

I 

9,  which  is  also  at  C^/U  *■  0.33.  The  loss  coefficient,  w,  is 

here  defined  as  the  (pitch  averaged)  total  pressure  difference  across 

blade  divided  by  the  local  inlet  dynamic  head.  For  the  solid  wall  the 


loss  profile  is  in  accord  with  what  one  would  expect,  with  low  losses 
near  the  midspan  region,  and  higher  values  at  the  two  endwall  regions, 
with  treatment,  however,  the  total  pressure  "loss"  coefficient  (as 
defined)  becomes  negative  from  midspan  to  hub  endwall,  showing 
clearly  that  the  hub  is  doing  work  on  the  flow. 

We  have  previously  discussed  the  differences  in  radial  streamline 
shifts  for  the  two  configurations.  The  radial  streamline  shifts 
derived  from  the  data  reduction  program  are  presented  in  Figure  10. 

In  the  figure  the  streamline  shift  between  the  RTE  and  STE  axial 
locations,  divided  by  span,  is  plotted  as  a  function  of  radius.  For 
the  untreated  case  there  is  a  large  streamline  shift  due  to  the 
blockage  occuring  near  the  hub,  while  for  the  treated  case  the  shift 
is  smaller  because  of  the  reduced  blockage. 

Using  this  information  about  the  radial  flows,  we  can  now  'Uuae 
the  loop  and  present  a  rotor  exit  total  pressure  profile  which  has 
been  convected  according  to  the  radial  streamline  shift  shown  in 
preceding  figure.  (Note  that  this  shifted  otal  pressure,  denoted  by 
RSH  was  the  reference  pressure  for  the  total  pressure  plots  presented 
in  Figure  6.)  Figure  11  shows  this  shifted  rotor  total  pressure 
profile  and  the  measured  midpitch  stator  total  pressure  profile,  for 
the  treated  hub  configuration,  at  Cx/U  =  0.33.  The  reference 
pressure  used  here  is  the  ambient  pressure,  and  the  dynamic  head  for 
nondimensionalization  is  based  on  the  midspan  stator  inlet  velocity. 

It  is  seen  that  in  the  neighborhood  of  the  midspan  the  two  profiles 
are  closely  the  same;  this  shows  the  consistency  between  the  data 
analysis  and  the  experimental  results  since  at  midpitch  the  flow  can 


be  considered  inviscid  and  the  total  pressure  should  be  constant 
along  a  streamline.  Note,  however,  that  near  the  hub  the  stator  total 
pressure  is  higher  than  the  shifted  rotor  total  pressure;  as  described 
previously  this  is  viewed  as  being  due  to  the  net  momentum  flux  out  of 
the  treatment,  which  is  seen  as  a  key  characteristic  of  hub  treatment. 

In  addition  to  the  data  presented  previously,  other  data  of  a 
diagnostic  nature  were  also  taken  to  investigate  the  sensitivity  of 
the  results  to  the  Reynolds  number.  Over  the  range  tested  (.75  x  10* 
to  1.35  x  10s)  there  was  little  change  (less  than  three  percent 
maximum)  in  any  of  the  measured  performance  parameters.  Since  the 
difference  between  treated  and  untreated  builds  at  any  of  the 
Reynolds  numbers  tested  was  an  order  of  magnitude  larger  than  this, 
Reynolds  number  effects  on  the  overall  phenomenon  are  small. 

Discussion 

From  the  experimental  results  it  appears  that  there  is  a  strong 
similarity  between  hub  treatement  and  casing  treatment;  flow  injection 
out  of  the  slots  and  consequent  energizing  of  the  endwall  region  and 
reduction  being  a  key  element  of  the  operation  of  both.  Accompanied 
with  the  casing  treatment  results,  it  is  also  believed  that  the 
hypothesis  of  "blade  stall"  and  "wall  stall"  can  be  a  useful  criterion 
for  successful  treatment  application.  From  a  more  detailed  fluid 
dynamic  point  of  view,  however  the  basic  mechanism  for  the  operation 
of  endwall  treatment  is  still  not  really  understood  in  a  quantitative 
manner.  In  this  regard,  the  terms  wall  stall  and  blade  stall  are  not 
really  descriptive  in  any  sort  of  fundamental  fluid  mechanic  sense. 


It  is  thus  important  to  quantify  more  precisely  the  features  of  the 
flow  regimes  that  are  associated  with  each  type  of  stall,  as  well  as 
the  detailed  effects  of  casing  treatment  on  the  stall  onset  process. 
The  present  set  of  experiments  can  therefore  be  viewed  more  as  an 
overall  examination  of  the  problem  and  it  is  clear  that  questions  such 
as  those  posed  in  ( 1 )  will  require  information  on  a  more  detailed 
level. 

There  is  one  further  point  that  can  be  mentioned.  As  previously 
discussed,  it  is  hypothesized  that  the  delay  in  stall  outset  due  to 
endwall  treatment  is  associated  with  momentum  transfer  between  slots 
and  blade  passage.  This  momentum  transfer  is  viewed  as  reducing  the 
endwall  boundary  layer  blockage.  We  have  explored  this  using  a  simple 
endwall  boundary  layer  flow  model  to  examine  the  effect  of  casing/hub 
treatment  on  boundary  layer  growth  (6).  It  is  found  that  the  momentum 
transfer  at  endwall  due  to  flow  in  and  out  of  the  grooves  has  a  much 
more  powerful  effect  on  the  boundary  layer  growth  than  the  wall  shear 
stress  has.  The  results  indicate  that  the  endwall  treatment  does  have 
a  very  strong  potential  to  decrease  the  boundary  layer  growth, 
although  the  calculated  effects  are  larger  than  those  observed  in  the 
limited  data  available. 
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Sumnary  and  Conclusions 

1.  An  investigation  has  been  carried  out  on  the  effects  of  a 
rotating  grooved  hub  treatment  on  stator  stall  and  performance.  Two 
different  sets  of  stators,  a  low  stagger  and  a  high  stagger  row,  were 
tested  with  and  without  the  axial  skewed  groove  hub  treatment. 

2.  It  was  found  that  for  the  low  stagger  row  the  application  of 
hub  treatment  did  not  have  a  major  effect  on  the  stator  performance. 
It  also  appeared  that  the  type  of  stall  encountered  with  this 
configuration  was  blade  stall. 

3.  For  the  high  stagger  row  the  application  of  axial  skewed 
slots  to  the  rotating  hub  under  a  compressor  stator  resulted  in  a 
substantial  change  in  both  the  stall  onset  point  and  peak  static 
pressure  rise  of  the  stator.  To  the  authors'  knowledge  this  is  the 
first  time  such  results  have  been  obtained. 

4.  It  is  hypothesized  that  the  success  of  the  treatment  was  due 
to  the  suppression  of  wall  stall  at  the  hub,  this  being  supported  by 
hot  wire  surveys  and  circumferential  traverses  at  stator  exit. 
Accompanied  with  the  results  of  the  low  stagger  case,  this  seems  to 
illustrate  again  that  endwall  treatment  will  be  effective  when  the 
type  of  stall  that  occurs  is  wall  stall. 

5.  It  is  seen  that  the  large  blockage  at  hub  endwall  region  is 
removed  when  the  treatment  is  applied.  The  high  losses  associated 


with  to  the  blockage  are  not  only  decreased,  but  the  total  pressure 
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actually  increases  across  the  stator  near  the  hub  region,  because  of 
the  work  done  by  the  (treated)  hub  on  the  main  flow. 
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Figure  3.  Stator  Static  Pressure  Rise  Characteristic  (low  stagger  stators) 
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Figure  7.  Stator  Inlet  and  Exit.  Flow  Angle  Distributions  with  Solid  Wall  (SW)  and  Hub  Treatment  (HT) 
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Figure  11.  Covected  Stator  Inlet  and  (Midpitch)  Stator  Exit.  Total  Pressure  Distributions 
(convection  of  inlet  profile  based  on  streamline  curvature  calculations). 


TASK  III:  INLET  VORTEX  FLOW  DISTORTION  IN  GAS  TURBINE  ENGINES 

During  the  past  year,  work  on  this  project  has  been  focussed  on 
the  completion  of  the  inlet  flow  facility  to  be  used  in  the  Wright 
Brothers  Wind  Tunnel,  the  basic  regimes  of  operation  of  this  facility, 
and  the  obtaining  of  initial  data  on  the  overall  flow  field  around  the 
inlet.  The  details  of  this  work  are  reported  in  GT  &  PDL  Report  No. 
166  (3.1),  which  is  essentially  the  M.S.  Thesis  of  Mr.  Wen  Liu. 

One  of  the  tasks  that  was  carried  out  in  connection  with  the 
facility  construction  was  the  fabrication  of  two  inlet  models,  one  of 
plexiglass  for  flow  visualization  and  one  of  aluminum  with  forty 
surface  pressure  taps.  The  inlets  were  designed  to  be  separation  free 
using  a  NASA  panel  method  code  and  semi-empirical  diffusion  factor 
separation  criterion  (3.3).  Another  task  was  the  construction  of  a 
shear  profile  generator,  which  creates  a  roughly  linear  velocity 
profile  over  the  central  fifty  percent  of  the  tunnel.  Also  included 
in  this  aspect  of  the  project  were  the  calibration  of  the  flow 
measuring  system,  the  survey  of  the  overall  tunnel  flow  quality,  and 
the  development  of  a  tripping  device  (wire)  on  the  inlet  to  ensure 
that  operation  was  not  in  the  subcritical  regime.  A  schematic  of  the 
inlet  model  is  shown  in  Fig.  3.1,  which  also  indicates  the  pressure 
tap  locations. 

Using  the  facility,  basic  parametric  information  has  been 
obtained  on  the  regimes  in  which  a  vortex  will  form.  This  information 
is  shown  in  Figure  3.2,  which  relates  to  an  inlet  at  ninety  degrees  of 
yaw,  in  a  flow  with  a  freestream  that  is  irrotational  upstream  of  the 
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inlet*  The  horizontal  axis  is  non-dimensional  inlet  height  above  the 
ground,  H/D,  where  D  is  the  inlet  inner  diameter,  and  the  vertical 
axis  is  the  velocity  ratio  It  can  be  seen  that,  as  has 

been  found  elsewhere:  (e.g.  3.3),  moving  to  larrer  values  of  H/D  or 
smaller  values  of  V^nie^/V®  decreases  the  tendency  for  an  inlet 
vortex  to  form.  This  type  of  parametric  data  has  also  been  obtained 
for  an  inlet  at  zero  degrees  of  yaw  in  a  flow  with  upstream  vertical 
vorticity,  as  shown  in  (3.1). 

Of  more  interest,  however,  is  the  data  on  the  static  pressure 
distribution  on  the  inlet  surface  for  an  inlet  at  ninety  degrees  of 
yaw.  These  data  are  plotted  in  terms  of  the  static  pressure 
coefficient,  Cp  (  =  [(Pi  -  Pref ) )/(?i  -  P)refl>,  versus 
azimuthal  position  round  the  inlet,  for  different  axial  locations,  and 
are  shown  in  Figure  3.3(a  to  h). 

The  Cp  distribution  is  shown  for  the  conditions  at  which  there 
is  a  strong  vortex;  H/D  «  1.0,  Vi/V»  »  16.5.  The  Cp  distribution 
versus  circumferential  location  (0  degrees  being  the  leading  edge,  and 
the  angles  varying  in  a  clockwise  manner  so  that  90  degree  is  top  and 
270  degree  is  bottom)  for  the  eight  axial  locations  along  the  inlet  is 
shown  in  Figures  3.3  (a)  to  3.3  (h).  Note  each  graph  is  plotted  on  a 
different  vertical  scale.  Figure  3.4  however,  is  a  composite  plot 
of  the  Cp  distributions,  to  the  same  scale. 

At  distances  far  away  from  the  inlet  lip,  the  Cp  distribution 
is  somewhat  similar  to  the  Cp  for  a  two  dimensional  cylinder, 
indicating  a  wake  region  behind  the  cylinder.  As  we  approach  the  lip, 
the  minimum  value  of  Cp  decreases  indicating  an  increase  in  flow 
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velocity.  The  shape  of  the  Cp  distribution  changes  gradually,  with 
a  region  of  high  static  pressure  appearing  at  the  lower  rear  of  the 
cylinder.  This  region  of  high  pressure  moves  to  about  the  260  degree 
position  at  axial  location  4  (1.0  D  from  the  lip).  (Note  that  there 
is  not  really  a  true  stagnation  point  due  to  the  component  of  surface 
flow  axially  along  the  inlet.) 

The  shift  in  the  location  of  the  low  velocity  region  from  the 
leading  edge  to  the  bottom  of  the  cylinder  suggests  a  change  in 
circulation,  increasing  towards  the  lip,  as  depicted  in  Figure  3.5. 

We  can  draw  an  analogy  to  the  study  reported  in  Goldstein  (3.4)  in 
which  a  cylinder  was  spun  to  create  circulation.  In  that  case  the 
front  stagnation  point  also  moved  towards  the  lower  part  of  the 
cylinder  as  circulation  increases. 

For  axial  locations  closer  to  the  lip,  the  Cp  plots  show  a  peak 
at  the  210  degree  location  indicating  a  region  of  low  velocity  there. 

A  tentative  explanation  for  this  is  as  follows.  From  the  flow  visual¬ 
ization  study  of  the  trailing  vortex,  it  was  shown  that  the  center  of 
the  vortex  trails  from  the  inlet  at  roughly  the  180  degree  position. 
Referring  to  Figure  3.6  we  can  see  that  there  is  an  interaction 
between  the  trailing  vortex  and  the  overall  flow  into  the  inlet,  which 
could  give  rise  to  a  region  of  low  velocity  at  about  the  210  degree 
location. 

In  addition  to  the  flow  associated  with  the  trailing  vortex,  the 
ground  vortex  which  is  below  and  in  front  of  the  lip  at  the  180  degree 
location  has  sane  influence  on  the  flow  into  the  inlet,  and  this  would 
also  be  expected  to  give  rise  to  a  low  velocity  region  near  the  210 
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degree  location,  and  a  high  velocity,  low  pressure,  region  near  the 
310  degree  location.  Because  the  ground  vortex  is  concentrated  in  a 
smaller  region  than  the  trailing  vortex,  the  flow  induced  by  the 
ground  vortex  is  locally  stronger  as  indicated  by  the  Cp  distribu¬ 
tions  at  axial  locations  7  and  8  (0.27  D  and  0.17  D  respectively  from 
the  inlet  lip);  the  velocity  at  the  310  degree  location  is  lower  than 
that  induced  by  the  trailing  vortex  at  about  the  150  degree  location. 
The  pressure  distribution  thus  appears  to  be  at  least  qualitatively 
explained  using  this  simple  model  of  varying  circulation. 

Ref  3.1  also  gives  details  of  the  flow  round  the  inlet  for  a  low 
value  of  where  we  obtain  a  totally  different  Cp 

distribution  near  the  inlet  lip.  The  composite  plot  (of  all  eight 
axial  locations)  for  this  situation  is  shown  in  Figure  3.7  where  we 
see,  as  the  inlet  lip  is  approached,  a  peak  in  static  pressure,  i.e.  a 
low  velocity  region  near  180  degrees.  This  may  be  explained  with 
reference  to  Figure  3.8,  which  we  regard  as  showing,  based  on  our 
previous  work  (3.5),  two  counter-rotating  tracking  vortices.  It  can 
be  seen  that  these  will  create  a  low  velocity  region  between  them, 
i.e.  near  180  degrees. 

In  summary,  in  both  situations  (the  high  Viniet  /V»  and  the 
low  V^nlet/»)  the  pressure  distributions  appear  to  be  in  accord 
with  the  conceptual  model  of  the  flow  round  the  inlet  that  we  have 


presented  in  our  previous  work 
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Figure  3.1.  Inlet  model  (D  =  0.15m) 
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Fig.  3.2  —  Regimes  of  vortex  formation  (inlet  at  90  degrees  of  yaw 
in  irrotational  freestream) 
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Fig.  3.4  Axial  Comparison  of  Circumferential  Cp  Distribution  of  Inlet  at 

90°  Yaw  Configuration  in  Irrotational  Upstream  Flow  at  High  V. 
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TASK  IV:  INVESTIGATION  OF  DISCRETE-BLADE  AND  3-D  FLOWS 
IN  HIGHLY-LOADED  TURBOMACHINES 

1 .  Blade  to  Blade  Effects 

The  following  Is  a  report  of  progress  made  in  the  past  year  In 
studying  blade-to- blade  (including  3-D)  effects  in  practical  turbo¬ 
machine  flows  of  high  loading.  Part  of  the  emphasis  is  on  contri¬ 
bution  to  improved  design  techniques ,  using  modified  non-linear 
analytic  techniques. 

Successful  non-linear  calculations  of  incompressible  three- 
dimensional  flow  through  high-swirl  blade  rows,  using  "very  thin" 
blades  have  been  reported  earlier.  Those  reults  confirm  expectations 
of  significant  3-D  effects.  But  they  also  suggested  the  importance  of 
exploring  additional  modifications  that  may  arise  from  blade-row 
blockage  ( blade  thickness ) . 

Proper  treatment  of  "leading-edge"  stagnation  points,  "Venturi" 
effects  in  the  intra-blade  region,  and  the  trailing  edge  Kutta 
conditions  could  all  be  expected,  it  was  thought,  to  modify  the  blade 
shapes  needed  to  accomplish  a  given  task.  Significant  differences  in 
required  camber,  depending  on  the  blockage,  was  to  be  expected. 

During  this  reporting  period  application  of  modified  non-linear 
wing  theory  methods  have  been  successfully  applied  to  treat  this 
problem. 

Taking  advantage  of  the  periodic  symmetry  of  typical  blade  rows, 
and  the  fact  that  they  frequently  comprise  many  blades,  an  asymptotic 
expansion  of  the  blade-row  "Green's  functions"  can  be  found.  Using 
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thls  result  in  a  more-or-less  classical  manner  for  the  solution  of 
fluid  flow  boundary-value  problems  leads  to  an  analytical  method  which 
we  call  "smoothing."  This  appears  to  greatly  reduce  the  work  required 
to  solve  blade  row  problems. 

This  "smoothing"  technique  has  now  been  successfully  applied  to 
the  "blockage"  or  "thickness"  problem  in  designing  highly-loaded 
blades  in  a  2-D  cascade  of  high  swirl. 

One  feature  of  the  above  results  in  that  one  can  now  explore, 
cheaply  and  efficiently,  a  wide  variety  of  intrablade  design  choices. 
Included,  for  example,  is  the  possibility  of  finding  out  what  blockage 
distributions  (chordwise)  are  favorable  or  unfavorable  for  given 
loading  distributions  (or  whirl  schedules).  This  work  comprised  the 
master's  thesis  of  Mr.  Thong  Dang. 
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2.  Non-Axisymmetric  Compressible  Swirling  Flow  in  Turbomachine  Annuli 

An  analysis  of  the  behavior  of  non-ax isymmetric  disturbances  in  a 
swirling  compressible  flow  typical  of  those  found  in  turbomachine 
annuli  has  been  completed.  The  results  of  this  work  will  appear  in 
the  near  future  as  GT  &  PDL  Report  No.  169. 

The  analysis  used  is  based  on  the  Clebsch  transformation.  For 
the  flows  analyzed,  the  different  types  of  flow  disturbances 
(vorticity,  entropy  and  pressure)  are  strongly  coupled,  due  to  the 
swirl.  Numerical  examples  are  presented  to  show  that  the  mean  flow 
Mach  number  can  have  an  influence  on  the  magnitude  and  downstream 
evolution  of  the  non-axisymmetric  disturbances.  In  particular,  at 
high  subsonic  Mach  number,  the  amplitude  of  the  disturbances  can 
oscillate  with  decreasing  modulation  in  the  flow  direction.  In 
addition  to  the  general  analysis  that  is  presented,  another  much 
simpler  model  (which  only  considers  one  class  of  non-axisymmetric 
disturbance)  is  also  used  to  examine  the  static  pressure  non¬ 
uniformities  that  can  occur  in  such  flows. 

To  illustrate  the  magnitude  of  the  effects  that  are  due  to 
compressibility,  fig.  4.1  shows  the  evolution  of  the  downstream  static 
pressure  distributions  for  three  different  far  upstream  Mach  numbers: 
0.1,  0.5,  and  0.9.  Note  that  the  example  shown  is  a  rotor  in  a 
annulus  of  constant  hub  and  tip  diameter.  With  the  present 
calculations,  it  is  necessary  to  set  the  far  upstream  Mach  number  at 
the  quite  high  value  of  0.9  to  obtain  a  downstreams  mach  number  in  the 
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high  subsonic  regime.  This,  however,  should  not  be  important  as  far 
as  the  downstream  evolution  is  concerned. 

3.  Theory  of  Blade  Design  for  Large  Deflections 

A  three  dimensional  method  of  designing  highly  loaded  blades 
(both  two-dimensional  and  annular  blades)  to  give  a  specified  distri¬ 
bution  of  swirl  has  been  formulated.  For  the  initial  work,  the  flow 
is  assumed  to  be  incompressible  and  inviscid  (the  annulus  has  constant 
hub  and  tip  radii),  the  blades  are  of  negligible  thickness;  and  a 
simple  free  vortex  swirl  schedule  is  assumed  for  the  case  of  annular 
blades,  numerical  calculations  have  been  carried  out  to  show  the 
magnitude  of  the  effects  in  parameter  regimes  of  interest.  The 
results  of  this  work  (two  papers)  have  been  submitted  to  the  ASME  for 
the  1983  Gas  Turbine  Conference.  Abstracts  of  the  two  papers  are 
reproduced  below.  The  papers  themselves  are  given  as  Appendix  B.  In 
addition,  the  theoretical  and  computational  details  will  appear  in  GTL 
Report  Mo. 168. 


Theory  of  Blade  Design  for  Large  Deflections 
Part  I:  Two  Dimensional  Cascades 

A  method  of  designing  two  dimensional  cascades  of  highly  loaded 
blades  to  give  a  specified  distribution  of  tangential  velocity  is 
presented.  The  fluid  is  assumed  to  be  incosg»ressible  and  inviscid  and 
the  blades  of  negligible  thickness.  The  blades  are  represented  by  a 
distributed  bound  vorticity  whose  strength  is  determined  by  the 
prescribed  tangential  velocity.  The  velocity  induced  by  the  bound 
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vortices  is  obtained  by  a  conventional  Biot-Savart  method  assuming  a 
first  approximation  to  the  blade  profile.  Using  the  blade  surface 
boundary  condition,  the  profile  is  then  obtained  by  iteration.  It  is 
shown  that  this  procedure  is  successful  for  large  pitch-chord  ratios 
and  large  deflections.  In  order  to  develop  a  method  for  use  in  three 
dimensions,  the  velocity  is  divided  into  a  pitchwise  mean  value  and  a 
value  varying  periodically  in  the  pitchwise  direction.  By  using 
generalized  functions  to  represent  the  bound  vorticity  and  a  Clebsch 
formulation  for  the  periodic  velocity,  series  expressions  are  obtained 
which  can  be  adapted  to  three  dimensional  problems.  Several  numerical 
results  are  obtained  using  both  approaches. 

Theory  of  Blade  Design  for  Large  Deflection 
Part  II:  Annular  Cascades 

A  method  of  designing  highly  loaded  blades  to  give  a  specified 
distribution  of  swirl  is  presented.  The  method  is  based  on  a  newly- 
developed  three-dimensional  analysis.  In  the  present  application,  the 
flow  is  assumed  to  be  incompressible  and  inviscid  (the  annulus  has 
constant  hub  and  tip  radii)  and  the  blades  are  of  negligible 
thickness.  A  simple  free  vortex  swirl  schedule  is  assumed.  The  flow 
velocity  is  divided  into  circumferentially  averaged  and  periodic 
terms.  The  Clebsch  formulation  for  the  periodic  velocities  is  used 


and  the  singularities  are  represented  by  periodic  generalized 
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functions  so  that  solutions  may  be  obtained  in  terms  of  eigen¬ 
functions.  The  blade  profile  is  determined  iteratively  from  the  blade 
boundary  condition. 

Results  from  the  computer  program  show  how  blade  number,  aspect 
and  hub-tip  ratios  affect  the  blade  shape.  The  blade  profiles  for  a 
given  swirl  schedule  depend  not  only  on  the  aspect  ratio  but  also  on 
the  stacking  position  (i.e.,  the  chordwise  location  at  which  this  thin 
blade  profile  is  radial),  and  so  too  do  the  mean  axial  and  radial 
velocities.  These  effects  occur  whether  the  number  of  blades  is  large 
or  small,  and  we  conclude  that  even  in  incompresible  flow  the  blade 
element  or  strip  theory  is  not  generally  satisfactory  for  the  design 
of  high  deflection  blades.  The  analysis  derives  conditions  for  the 
blade  profiles  on  the  walls  of  the  annulus  as  well  as  the  mean  swirl 
velocity  there  and  at  the  trailing  edge. 

In  the  limit  when  the  number  of  blades  approaches  infinity,  a 
bladed  actuator  duct  solution  is  obtained.  The  conditions  for  the  blade 
profile  at  the  walls  are  absent,  but  the  stacking  position  and  aspect 
ratio  still  affect  the  axial  and  radial  velocity  distributions  for  the 


8 ante  swirl  schedule 


4 


Numerical  Studies  of  Secondary  Flow  in  a  Bend  Using  Spectral 
Methods 


A  numerical  algorithm  for  solution  of  the  incompressible  Navier- 
Stokes  equations  for  flow  in  a  bend  of  square  cross-section  has  been 
formulated.  The  momentum  equation  has  been  written  in  rotational  form 
as 

|2«vxu-Vs+  vV*V  (1) 

3t 

and  the  continuity  equation  is  simply 

V  •  T  -  0  (2) 

For  Eq.  1  and  2,  is  the  velocity  vector,  3  the  vorticity,  tt 
the  stagnation  pressure  normalized  by  the  density,  and  v  the  kinematic 
viscosity.  A  cylindrical  coordinate  system  (r,  B,  z)  is  used  with  the 
square  section  of  the  bend  lying  on  the  r  -  z  plane. 

The  fractional  time  stepping  technique  has  been  proposed  to  advance 
the  Navier-Stokes  equation1  in  time  as  follows: 

(1)  The  convective  (or  advective)  part: 


The  first  fractional  step  is  the  convective  part.  A  semi-implicit 
Crank-Nicholson-Adams-Bashforth  Scheme  is  used  and  yields 
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vg  is  the  value  of  the  bend  through-flow  velocity  at  the 
inflow  boundary.  Superscript  (1)  denotes  the  first  fractional 
time  step  while  superscript  n  denotes  the  time  level.  For  this 
part,  no  outflow  B.C.  are  required  as  the  spectral  method  does 
not  introduce  spurious  B.C.  The  inflow  boundary  conditions  are: 


vr  «  0 

v q  -  CQ(1  -  R2) (1  -  Z2) 

V2  -  o 

(2)  The  pressure  correction  step: 


(4) 


The  second  fractional  step  enforces  the  incompressibility  condition 
and  the  normal  wall  boundary  conditions,  which  are  not  required  by 
the  first  fractional  step.  This  is  achieved  via  the  solution  of  v 
from 


L  [w] 


l_  7  .  ; 

At 


(5) 


where  L  is  the  appropriate  Laplacian  operator  in  (R,$ ,Z)  space 
mapped  from  the  cylindrical  co-ordinate  space  (r,  0,  z).  The 
boundary  conditions  imposed  on  *  are: 
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2  At  V  2  2 
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— 
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at  the  outflow  boundary 

(6d) 

3  ♦ 

(8o  Is  the  angle  o€  the  bend,  rl  and  ro  are  the  inner  and 
outer  radius  of  the  bend) 

(Note  that  the  imposed  outflow  boundary  condition.  Equation  (6d), 
is  not  strictly  correct,  and  we  will  be  attempting  to  improve  (6d) 
utilizing  a  linearized  secondary  flow  analysis.)  in  eqn.  (6),  we 
have  used  subscript  (r,9,z)  to  denote  vector  component  in  the  r- , 
9-,  z-  direction.  The  corrected  velocity  field  g(2)n+1  j_8 

obtained  through 

$(2)n+1  .  $<1)n+1  _  4t7¥  (7, 

(3)  The  viscous  correction  stepi 

The  final  fractional  step  consists  of  ingxssing  the  no-slip 
conditions  at  the  walls  together  with  appropriate  inflow  and 
outflow  B.C.  The  beclcward  Euler-like  aches*  is  used  and  it  yields 
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Steps  (1)  and  (2)  are  of  0(At2)  while  step  (3)  is  of 
O(vAt).  However,  due  to  the  curvature  term,  the  pressure  and  the 
viscous  operator  do  not  commute,  therefore  the  overall  accuracy 
will  be  of  0 (At) •  However,  0(At2)  nay  be  achieveable  by 
alternating  the  direction  of  fractional  stepping.  It  should  also 
be  mentioned  that  the  implementation  of  fractional  step  ( 1 )  and 
(2)  only,  will  yield  an  inviscid  solution  of  the  flow. 

i 

t 

We  have  proposed  to  expand  the  velocity  field  $  and  the 
stagnation  pressure  w  in  triple  chebyshev  series.  Consequently,  we 
have  developed  software :  ( 1 )  to  compute  and  invert  single ,  double  and 

triple  Chebyshev  series  efficiently  using  FFT;  (2)  to  compute  the 
derivatives  of  function  efficiently  given  the  Chebyshev  coefficients 
of  the  function.  This  includes  one-dimensional,  two-dimensional  and 
three-dimensional  functions.  (3)  To  solve  the  three-dimensional 
Poisson  equation  in  a  cylindrical  co-ordinate  system  with  both  Neumann 
type  of  boundary  conditions  (pressure  step),  and  Oirichlet  boundary 

i 

conditions  (various  step).  We  have  tested  the  software  developed  in 

t 

I 

i 

! 

i 


(1)  and  (2)  rather  exhaustively  (against  known  analytical  functions) 
to  ensure  the  correctness  of  the  software. 

For  the  three-dimensional  Poisson  solver,  we  have  initially 
proposed  the  use  of  a  tensor-product  technique  with  diagonalization  in 
the  4  and  z  direction  and  spectral  iteration  in  the  radial  direction. 
The  spectral  iteration  operator  chosen  is  the  second  finite  difference 
operator  for  3* /Sr*  +  (1/r)  3/3r  with  the  appropriate  boundary 
conditions.  However,  we  are  now  in  favor  of  a  slightly  different 
method  for  the  solution  of  three-dimensional  Poisson  equation  in  the 
cylindrical  co-ordinate  system.  In  this  method,  we  still  perform  the 
diagonalization  in  the  ♦  -  and  z-  direction,  but  we  now  use  a  collo¬ 
cation  approximation  for  the  operator  in  the  radial  direction.  An 
eigensolver  is  then  used  to  diagonalize  the  collocation  operator  to 
reduce  the  storage  and  number  of  operations  required  for  the  solution. 
The  current  technique  still  requires  a  preprocessing  stage.  The 
software  has  been  developed  for  this  and  has  been  tested  on  model 
equations  with  homogenous  B.C.  with  known  analytical  solution.  For 
inhomogenous  boundary  conditions,  it  is  necessary  to  invent  a  function 
(usually  a  polynomial)  to  make  the  boundary  conditions  homogenous  with 
a  resulting  modification  of  the  source  term  on  the  R.H.S.  of  the 
Poisson  equation  for  this  method  to  work.  We  have  found  that  the 
invention  of  such  a  function  to  be  rather  inconvenient  (especially  for 
the  case  of  inhomogenous  Neumann  Boundary  condition  which  is  the  case 


for  the  Pressure  step  here).  We  expect  to  modify  our  current  scheme 


so  that  the  case  of  a  general  boundary  condition  of  the  mixed  type  can 
be  handled  conveniently. 

It  is  noted  that  the  use  of  Chebyshev  polynomials  in  the  r-z 
plane  leads  to  a  good  resolution  of  the  boundary  layers  at  the  walls. 
However,  its  use  in  the  $  direction  does  lead  to  some  wasteful  reso¬ 
lution  in  the  inflow  and  outflow  boundaries.  Besides,  since  spectral 
methods  are  global,  therefore  any  error  incurred  during  the  imposition 
of  outflow  boundary  conditions  (e.g.  6d)  may  not  tend  to  be 
localized.  Thus  we  may  anticipate  the  modification  to  finite 
difference  approximation  in  the  4 -direction  to  alleviate  the  wasteful 
resolution  at  the  inflow  and  outflow  boundaries  and  to  localize  the 
errors  incurred  in  the  imposition  of  outflow  boundary  condition  such 
as  that  in  equation  6d. 

5.  Numerical  calculation  of  flow  in  casing  treatment  grooves 

Work  is  in  progress  on  the  calculation  of  the  flow  in  a 
rectangular  slot,  characteristic  of  those  used  in  rotor  casing 
treatment.  The  flow  is  assumed  to  be  inviscid,  incompressible, 
irrotational ,  and  two-dimensional.  The  boundary  condition,  in 
addition  to  the  usual  vanishing  of  the  normal  velocity  at  the  solid 
walls,  is  the  time-dependent  pressure  distribution  imposed  by  the 
passage  of  rotor  blades  on  the  open  end  of  the  rectangular  slot.  In 
mathematical  terms,  we  solve  the  following  problem: 
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in  -1  <  x  <  1  ,  and  -1  <  y  <  1  (see  fig.  4.2) 

subjected  to  the  boundary  conditions: 


34 

= 

0 

at  x  =  ±  1 

(2a) 

3x 

34 

= 

0 

at  y  =  -  1 

(2a) 

3y 

34 

2 

and  ft 

It 

+  (V4) 

+ 

p  =  constant  at  y  =  1 

(2c) 

The  above  equations  have  been  appropriately  made  dimensionless: 

4  is  the  velocity  potential  and  p  is  the  pressure  distribution  at  the 
open  end.  Q  (  2  x  blade  passing  frequency  x  slot  length/blade  tip 
speed)  can  be  taken  as  a  reduced  frequency  parameter  providing  a 
measure  of  the  importance  of  unsteadiness. 

A  typical  time-dependent  pressure  distribution  at  8  different 
time  levels  at  the  open  end  of  the  groove  constructed  from  the  blade 
pressure  distributions  given  in  the  Ph.D.  thesis  by  Smith  (4.2),  is 
shown  in  fig.  (4.3a).  The  time-averaged  pressure  distribution  is 
shown  in  fig.  (4.3b). 

Several  numerical  methods,  including  spectral  and  finite 
difference  schemes,  were  tried.  Spectral  expansions  were  either  in 
Chebyshev  polynomials  or  cosine  functions  (the  natural  eigenfunction 
of  the  system) .  Temporal  resolution  was  either  in  Fourier  series 
(taking  advantage  of  the  time  periodicity)  or  in  simple  finite 


difference  form  (Euler,  leapfrog,  etc.).  Because  of  the  nonlinearity 
of  the  boundary  condition  at  y  ■  1  (eg.  2c),  iteration  of  the 
relaxation  type 

-  4*  +  U  [a  04/31)  +  (R4)2  +  p  -  constant]  (3 

was  used  in  the  spectral  calculation.  Satisfactory  convergence  was 
only  obtainable  for  imaginary  relaxation  factor  u  and  5  spectral 
modes,  neglecting  the  D.C.  (timewise)  component  of  the  pressure  at 
the  open  end.  The  use  of  spectral  modes  greater  than  5  yielded  a 
n on-satisfactory  convergent  solution.  For  real  relaxation  factor  u , 
the  computational  procedure  was  unstable.  For  the  finite  difference 
schemes  that  were  tried,  no  satisfactory  convergent  solution  is 
obtained.  The  details  and  results  of  work  described  above  is 
documented  in  the  M.s.  Thesis  of  P.  Kletzkine. 

Recently,  we  have  implemented  the  calculation  of  the  time- 
averaged  flow  using  the  pressure  distribution  shown  in  fig.  4.3b.  The 
spectral  expansion  used  for  4  is  a  cosine  series.  The  nonlinearity  of 
the  boundary  condition  at  the  open  end  (fig.  4.3c)  results  in  multiple 
solutions  for  the  same  given  pressure  distribution.  The  most 
physicaliy  realistic  solution,  corresponding  to  the  pressure 
distribution  in  fig.  4.3b,  is  shown  in  fig.  4.4. 

We  are  currently  attempting  to  develop  a  roburst  technique  for 
the  solution  of  equation  1  subjected  to  boundary  conditions  2. 
Currently,  we  are  also  investigating  other  ways  (e.g.  using  the 
primitive  variables  p  and  velocity  components  u  and  v)  of 


numerically  computing  such  a  flow  field.  For  this  latter  work,  the 
flow  does  not  have  to  be  irrotational . 


4.1  Kells, 
Coette 
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GENERAL  PROGRESS  CM  WRAPT 


During  the  last  year  the  Air  Force  OSR  sponsored  educational 
program,  called  Air  Force  Research  in  Aeronautical  Propulsion 
Technology  (WRAPT),  was  initiated.  An  industry-faculty  steering 
committee  was  established.  It  consists  of 


Dr.  Michael  Salkind 


AFOSR,  Chairman 


Mr.  Walter  Schrader 


AVCO  Lycoming 


Dr.  Lynn  Snyder 


Mr.  Elmer  Wheeler 


Garrett  Aerospace  Corp. 


Dr.  Fred  Erich 


GE  Lynn 


Dr.  David  Wisler 


GE  Evendale 


Mr.  Morris  Zipkin 


Pratt  and  Whitney  (GPD) 


Professor  Sanford  Fleeter  Purdue  University 


Professor  Peter  Jenkins 


Texas  ASM  University 


professor  E.  E.  Covert 


Massachusetts  Institute  of 
Technology,  Secretary 


This  group  has  met  three  times,  established  procedures,  put  out 
the  first  year's  brochure  and  generally  put  the  program  in  motion. 

This  spring  22  applications  were  received,  9  offers  were  made 
and  5  acceptances  have  been  received. 

The  Steering  Committee  met  to  review  the  program  for  next  year 


and  decided  to  continue  the  program  in  the  present  form. 


< 
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in  turbomachines,  etc.),  presented  at  Nanjing  Aeronautical 
Institute,  Nanjing,  China,  and  Beijing  Institute  of  Aeronautics 
and  Astronautics,  Beijing,  China,  October  1981. 

E.M.  Greitzer,  seminar  at  NASA  Lewis  Research  Center 
"Non-Recoverable  Stall  in  Axial  Compressors,"  January,  1982. 

E.M.  Greitzer,  seminar  at  NASA  Lewis  Research  Center, 
"Mechanisms  of  Inlet  Vortex  Formation,"  June,  1982. 

E.M.  Greitzer,  lecturer  in  Southwest  Applied  Mechanics  Lecture 
Series  "The  Stability  of  Pumping  Systems",  October  1982. 

A.H.  Epstein,  "Time  Resolved  Measurements  in  a  Transonic 
Compressor, "  paper  presented  at  1982  ASME  Gas  Turbine 
Conference,  London,  England,  April,  1982. 


As  described  in  previous  reports,  there  are  strong  interactions 
between  the  Gas  Turbine  Laboratory  and  the  aircraft  engine  industry, 
as  well  as  Air  Force  and  other  government  laboratories.  In  addition 
we  have  an  active  seminar  program  to  bring  speakers  from  industry 
and/or  government  to  M.l.T.  During  the  past  year  these  have  included 


Dr.  J.  Caspar,  United  Technologies  Research  Center 
"Numerical  Calculation  of  Cascade  Flows  Using  Non-Orthogonal 
Grids" 

Mr.  R.B.  Davis,  Pratt  and  Whitney  Aircraft,  Government  Products 
Division 

"Designing  for  Transients  in  High  Performance  Gas  Turbine 
Engines" 

Dr.  M.E.  Goldstein,  NASA  Lewis  Research  Center 
"Generation  of  Tollmien-Schl ichting  Waves  by  Free-stream 
Disturbances" 

Dr.  ft.  Mani,  Corporate  Research  Center,  General  Electric  Company 
"Some  Aspects  of  Gas  Turbine  Engine  Post-Stall  Behavior" 
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Dr.  R.  Norton,  Rolls  Royce  U.S.A.  (on  leave  from  M.l.T. 
Department  of  Aero/ Astro) 

"Transonic  Flow  in  a  Nozzle  Guide  Vane" 

Mr.  R.  Novak,  Aircraft  Engine  Group,  General  Electric  Company 
"A  Mixed-Flow  Cascade  Passage  Design  Procedure  Based  on  A  Power 
»  Series  Expansion" 

Mr.  B.A.  Robideau,  Pratt  &  Whitney  Aircraft,  Commercial  Products 
3  Division 

"Practical  Application  of  Aerodynamic  Principles  to  Multistage 
Compressor  Design" 

Dr.  M.  Sajben,  McDonnell-Douglas  Corporation 

"Unsteady  Aspects  of  Internal  Flows  with  Shock/Boundary  Layer 
Interaction" 

Dr.  M.J.  Werle,  United  Technologies  Research  Center 
"Boundary  Layer  Separation  in  Turbomachinery" 
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6.  DISCOVERIES,  INVENTIONS  AND  SCIENTIFIC  APPLICATIONS 

Discoveries 

We  have  demonstrated  (for  the  first  time)  that  compressor  hub 
treatment  can  markedly  improve  the  stall  margin  of  an  axial 
compressor  stator  row. 

Inventions  and  Scientific  Applications 

A  transonic  flow  computation  procedure,  having  stagnation 
pressure  errors  on  the  order  of  1%,  has  been  developed  for 
turbomachinery  cascades. 
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7.  CONCLUSIONS 

In  brief,  although  some  of  the  projects  appear  farther  along  than 
others,  the  M.I.T.  Multi-Investigator  effort  appears  to  be  basically 
on  target.  In  addition,  it  now  appears  that  we  are  seeing  "build-up" 
efforts  of  the  previous  years  start  to  bear  fruit,  as  reflected  in  the 


publications  produced 


